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FOREWORD 

This Indian Standard was adopted by the Bureau of Indian Standards, after the draft finalized by the Ground 
Water and Related Investigations Sectional Committee had been approved by the Water Resources Division 
Council. 

Groundwater is available almost everywhere. However, its distribution is not uniform due to varying 
hydrogeological, topographical and climatic conditions. As a result, groundwater is not always available in the 
required quantity and/or quality, particularly in hard rock terrains where the fractures and weathered parts are the 
only conduits for groundwater. Therefore, collection of information on prospective groundwater zones, though 
being a costly affair, is an essential prerequisite. Surface geophysical methods are currently recognized as cost 
effective techniques that are useful for collecting this kind of information. Measuring physical properties of the 
earth and their variation, and then associating finally them with hydrogeological characteristics is the overall 
domain of groundwater geophysics. 

Of the various geophysical techniques available today, the electrical resistivity method is probably most commonly 
used due to its relatively simple and economical field operation, its effective response to ground water conditions, 
and the relative ease with which interpretations can be made. This type of survey is occasionally supplemented 
by other techniques such as induced polarization, spontaneous potential, and Mise a la Masse galvanic electrical 
techniques. Other geophysical methods in order of preference used for hydrogeological purpose are 
electromagnetic, refraction seismic, magnetic, gravity and seismic reflection surveys. More recently developed 
geophysical techniques include ground probing radar, electrokinetic sounding, and nuclear magnetic resonance, 
but these methods are not in widespread use and are not considered further in this report. 

Because surface geophysical surveys are carried out at the surface of the earth, the responses received from 
different depths often lack unique characteristics. That is, ambiguity exists in interpreted results and the effective 
application of these methods often depends on the skill and experience of the investigator, knowledge of the 
hydrogeological conditions, and the usefulness (and limitations) of the technique(s) themselves. The application 
of two or more geophysical techniques may also be a useful approach to use in some field surveys. Integration of 
information received from other scientific surveys, such as remote sensing, hydrogeologic characterization, 
chemical analysis of well water samples, etc, is also useful for interpreting the filed data. 

Modern geophysical techniques are highly advanced in terms of instrumentation, field data acquisition, and 
interpretation. Field data are digitized to enhance the signal-to-noise ratio, and computers are used to more 
accurately analyze and interpret the data. However, the present-day potential of geophysical techniques has 
probably not been fully realized, not only because such surveys can be expensive, but also because of the inadequate 
understanding of the application of relevant techniques in contrasting hydrogeological conditions. 

It has been assumed in the formulation of this standard that the execution of its provisions is entrusted to 
appropriately qualified and experienced people, for whose guidance it has been prepared. 

In reporting the results of a test or analysis made in accordance with this standard, if the final value observed or 
calculated, is to be rounded off, it shall be done in accordance with IS 2 : 1960 'Rules for rounding off numerical 
values (revised)' . 
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Indian Standard 

SURFACE GEOPHYSICAL SURVEYS FOR 
HYDRO GEOLOGICAL STUDIES 



1 SCOPE 

The application of surface geophysical methods is an 
evolving science that can address a variety of objectives 
in groundwater investigations. However, because the 
successful application of geophysical methods depends 
on the available technology, logistics, and expertise of 
the investigator, there can be no single set of field 
procedures or approaches prescribed for all cases. 
Accordingly, this standard described guidelines that 
should be useful for conducting geophysical surveys 
for a variety of objectives (including environmental 
aspects), within the limits of modern-day 
instrumentation and interpretive techniques. The more 
commonly used field techniques and practices are 
described, with an emphasis on electrical resistivity, 
electromagnetic, and seismic refraction techniques as 
these are widely used in groundwater exploration. 
Theoretical aspects and details of interpretational 
procedures are referred to only in a general way. 

2 REFERENCES 

The following standards contain provisions which 
through reference in this text, constitute provisions of 
this standard. At the time of publication, the editions 
indicated were valid. All standards are subject to 
revision, and parties to agreements based on this 
standard are encouraged to investigate the possibility 
of applying the most recent editions of the standards 
indicated below: 

IS No. Title 

15681 : 2006 Geological exploration by 
geophysical method (seismic 
vibration) — Code of practice 

15736 : 2006 Geological exploration by 
geophysical method (electrical 
resistivity) — Code of practice 

3 TERMINOLOGY 

For the purposes of this standard, the following terms 
and definitions shall apply. 

3.1 Acoustic Impedance — Product of seismic 
velocity and density of a layer. Reflection of seismic 
wave depends on contrast in acoustic impedance. 

3.2 Anisotropy — Variation in physical property with 
direction of measurement is anisotropy. In electrical 
resistivity method micro-, macro- and pseudo- 
anisotropy are involved. Anisotropy of a geoelectrical 



layer is given as A. = ^/p t /p L where p t and p L are 
transverse and longitudinal resistivities of a layer. 

3.3 Apparent Resistivity — Ratio of measured voltage 
to input current multiplied by geometric factor of 
electrode configuration. It would be true resistivity if 
the subsurface is homogeneous (scale of homogeneity 
referred to dimension of electrode geometry). 

3.4 Aquifer — Formation or group of formations or a 
part of formation that contains sufficient permeable 
material and is saturated to yield significant quantities 
of water to wells and springs. 

3.5 Blind Zone — Layer having seismic velocity less 
than that in the layer overlying it . 

3.6 Bouguer Correction — Correction made in 
observed gravity data to account for the attraction 
(gravitational) of the rock between the datum and the 
plane of measurement. It is 0.041 85 p h mgal, where p 
is the density of the rock between the datum and the 
plane of measurement and h is the difference in 
elevations between the datum and the plane. 

3.7 Bouguer Anomaly — Anomaly obtained after 
applying latitude, terrain, and elevation (free air and 
Bouguer) corrections to the observed gravity value and 
finally subtracting it from measured value at some 
particular station in the survey area. 

3.8 Contact Resistance — Electrical resistance 
developed between an electrode planted in the ground 
and the ground material immediately surrounding it. 
Contact resistance is reduced by putting water at the 
electrodes. 

3.9 Convolution — Defined as the integral of the product 
of the two functions after one is reversed and shifted. In 
digital signal processing, frequency filtering can be 
simplified by convolving two functions (data with a filter) 
in the time domain, which is analogous to multiplying 
the data with a filter in the frequency domain. 

3.10 Dar Zarrouk Parameters — Longitudinal unit 
conductance (S) and transverse unit resistance (I) of a 
geoelectrical layer. These are defined as 

S = I,h/p = (h l /p l + h 2 /p 2 +h 3 /p 3 + ...) 
T = T,hp = (h 1 p l + h 2 p 2 + h 3 p 3 + ...) 

where h u h^, h 3 are thickness and p 1; p 2 , p 3 .... are 

electrical resistivity of different subsurface layers. 
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3.11 Deconvolution — Process of inverse filtering to 
nullify the undesired effect of an earlier filter operation. 

3.12 Dipole-Dipole Electrode Configuration — 

Configuration in which the spacing between current 
electrode pair and that between potential electrode pair 
is considerably small in comparison to the distance 
between these two pairs. As the current and potential 
electrode pairs are moved further apart, the depth of 
exploration is progressively increased. The dipole- 
dipole configuration can be azimuthal, equatorial, 
radial, parallel, axial and perpendicular. Geometric 
factors are : 

2nr 3 /L I sin ... azimuthal, 2k r 3 ILl ... equatorial, 

n r'lL / cos 9 . . . radial, 2%r 3 /L I ( 3 cos 2 9 - 1 ) 

. . . parallel, 

K r'lL I ... axial and 2nr 3 /3L I ( sinB cosB) 

. . . perpendicular 

L is length of current dipole, / is length of potential 
dipole, r is distance between centres of current and 
potential dipoles, is the angle between the two dipole 
axes. 

3.13 Diurnal Correction — Correction applied to 
magnetic data to compensate for daily fluctuations of 
the geomagnetic field. 

3.14 Drift Correction — Quantitative adjustment to 
account for a uniform change in the reference value 
with time. Drift in gravity meters is mainly due to creep 
in the springs of the gravimeter. Correction to measured 
values are made by repeating readings at 3 h to 4 h at a 
fixed station. 

3.15 Eddy Current — Current induced in a conductive 
body by the primary electromagnetic (EM) field. The 
secondary EM field produced by the eddy current 
opposes the primary field. 

3.16 Equivalence — If target response is a function 
of product or ratio of two parameters (say bed thickness 
and resistivity), variation in the parameters keeping the 
ratio or product constant can yield almost same 
response and the various combination of parameters 
are said to be equivalent. This brings in ambiguity in 
parameter estimation. It is pronounced if the target is 
buried and relatively thin. In multi-layer geoelectrical 
sequence the intermediate layers show equivalence 
over a range of parameters. 

3.17 Free-Air Correction — Correction applied to 
gravity data to account for the fact that gravity 
measurements are made at different distances 
(elevations) from the center of the earth. The correction 
value is 0.308 6 h mgal, where h is the difference 
between the elevation of the datum and the plane of 
measurement. Free-air gravity anomaly is obtained 



after applying correction for the latitude and elevation. 

3.18 Geoelectrical Layer — Layer having 
characteristic of uniform electrical resistivity. 

3.19 Geometric Factor — Numerical value dependent 
upon the arrangement of electrodes which when 
multiplied by the measured voltage-to-current ratio 
gives the apparent resistivity. 

3.20 Geophone — Instrument which detects seismic 
energy and converts it into electrical voltage. Relative 
motion between a suspended coil and a magnet due to 
seismic wave generates a voltage in the coil whose 
amplitude is proportional to the velocity of the exciting 
seismic disturbance. 

3.21 Generalized Reciprocal Method — It is a 

technique wherein in-line seismic refraction data 
consisting of forward and reverse travel times are used 
for delineating undulated refractors at a depth. The 
travel times at two adjacent geophones are used in 
refractor velocity analysis and time-depth calculations. 
At the optimum inter-geophone spacing, the upward 
traveling segments of the rays to each geophone emerge 
from near the same point on the refractor. The depth 
conversion factor is relatively insensitive to dip angles 
up to about 20°, because both forward and reverse data 
are used. As a result, depth calculations to an undulating 
refractor are particularly convenient even when the 
overlying strata have velocity gradients. The GRM 
provides a means of recognizing and accommodating 
undetected layers, provided an optimum inter- 
geophone spacing value can be recovered from the 
travel-time data, the refractor velocity analysis, and/ 
or the time-depths. 

3.22 Gradient Configuration — A variation of the 
Schlumberger configuration where the current 
electrodes (AB) are kept at infinity, that is at a large 
separation and central l/3rd space is scanned by a small 
potential dipole (MN). The geometric factor is n/MN 
(AB/2) 2 { 1- x 2 /(AB/2) 2 } 2 I { l+x 2 l{ABI2) 2 } where x is 
the distance between the centre of the configuration 
and the centre of the potential dipole. 

3.23 Half-Schlumberger Configuration — 

Configuration in which one of the current electrodes 
is kept at infinity (large distance) and need not be 
collinear with the other three electrodes. It can be used 
for soundings along radial lines. The apparent 
resistivity is given as p a = 2na 2 /I (dV/da), where a is 
the distance between the active current electrode and 
center of the potential electrode spacing, da is the 
potential electrode spacing and dV is the potential 
difference. 

3.24 Homogeneity — Characteristic of a formation 
with uniform physical property or properties. It is a 
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function of scale of measurement in relation to the 
uniformity in physical property. Inhomogeneity or 
heterogeneity indicates non-uniformity or dissimilarity 
in physical property with reference to the scale of 
measurement. 

3.25 In-phase — Component of a secondary EM field 
with the same phase angle as that of the exciting 
primary EM field; that is, in-phase component attains 
maxima and minima in step with the primary field. 

3.26 Lee-partitioning Configuration — A variation 
of the Wenner array where one additional electrode is 
placed at the centre between the potential electrodes. 
Potential difference between the central electrode and 
either of the two other potential electrodes is measured. 
The geometric factor is 4na. Identical values of 
potential difference at each separation indicate that the 
sampled ground is homogeneous. 

3.27 Longitudinal Conductance — Ratio of the 
thickness of a geoelectric layer to its resistivity 
[conventionally expressed as S = hip (mhos)]. 

3.28 Magnetic Permeability — Ratio of magnetic 
induction (flux density) in a body to the strength of 
the inducing magnetic field. 

3.29 Magnetic Susceptibility — Ratio of the intensity 
of magnetization produced in a body to the strength of 
the magnetic field. 

3.30 Migration — That part of processing of seismic 
reflection data required to plot the dipping reflections 
at their correct position. 

3.31 Non-polarizing Electrode — Electrode which is 
not affected by electrochemical potential generated 
between the electrode and ground material in which it 
is planted. A copper rod placed in a copper sulphate 
solution contained in a porous ceramic pot is 
commonly used as a non-polarizing electrode. 

3.32 Normal Moveout — The effect of variation of 
shot-geophone distance on time of arrival of seismic 
reflection. 

3.33 Off-set Wenner Configuration — A modification 
in Wenner configuration to remove or minimize the 
effect of lateral inhomogeneities. In this configuration 
five equally spaced collinear electrodes are planted in 
the ground. Average is taken of consecutive normal 
Wenner measurements taking the left four and right 
four electrodes. 

3.34 Overburden — That part of the host medium 
which lies above the target and is usually of no interest 
in exploration, but has physical properties that affect 
the measurements. 

3.35 Phasor Diagram — Graph obtained by plotting 



in-phase and quadrature components of secondary EM 
field for different frequencies of primary field. The 
values of in-phase and quadrature components are 
plotted along x and y axes, respectively. Theoretical 
phasor diagrams are generated for different 
conductivity ratios and ratios of layer thickness to 
transmitter-receiver coil separation, and field data plot 
is matched. 

3.36 Plus-Minus (Hagedoorn) Method — Used to 
interpret seismic refraction data. The method uses 
reversed refraction profiles with shots at opposite ends 
and the addition and subtraction of travel times for 
various locations between the shots to give the depth 
to the refractor and its velocity. 

3.37 Polar Diagram — Method of plotting resistivity 
sounding data. The apparent resistivity values of the 
radial soundings conducted at a point are plotted for 
various current electrode separations. Results can be 
used to infer fracture orientations. 

3.38 Proton Precession Magnetometer — It is also 
known as nuclear precession magnetometer. Because 
of spin, proton has a magnetic moment. The axes of 
precession are oriented randomly. A magnetic field 
normal to the earth's magnetic field polarizes the nuclei 
for a short period and a voltage at precession frequency 
is induced in a measuring coil which indicates the value 
of earth's magnetic field at the point of measurement. 

3.39 Quadrature — Out-of-phase or imaginary 
component of secondary EM field, it is the component 
which is 90° out of phase with the inducing primary 
EM field. The ratio of the strengths of in-phase and 
quadrature components of secondary EM fields 
indicate conductivity characteristics of the target. 

3.40 Reflector — Interface which separates two layers 
of contrasting acoustic impedance giving rise to 
reflection. 

3.41 Refractor — Layer along which the refracted 
(head wave) wave travels at a velocity that is higher 
than that in the overlying layer. 

3.42 Remanent Magnetization — In-situ residual 
magnetization remaining in rock after removal of 
inducing field. 

3.43 Schlumberger Configuration — Collinear four- 
electrode configuration of current and potential 
electrodes in which potential electrodes are kept close 
to the center of the configuration. Conventionally, the 
separation between potential electrode (MAO is less than 
1/5 of the current electrode separation (AB). The 
geometric factor is Jt {{ABI2) 2 -(MNI2) 2 }IMN. 

3.44 Skin Depth — Effective depth of penetration of 
EM field in a medium. Skin depth is defined as the depth 
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where EM field intensity reduces to about 37 percent of 
its original value at the surface of the earth. It is dependent 
upon the conductivity and magnetic permeability of the 
medium, and the frequency of the EM field. Increases 
in these values reduce the skin depth, as does the 
presence of conductive overburden. It is expressed as 

z = 500 J~ (p//), where p is resistivity of the ground and 

/ is the frequency of the EM field generated. 

3.45 Snell's Law — When a seismic wave is incident 
at a particular angle (/) on a boundary between two 
media (having different seismic velocities v i and v 2 
v 2 > v,), the wave gets refracted at an angle (r) at the 
boundary according to Snell's law which states that 
sin ilv i = sin r /v 2 . 

3.46 Stacking — Process of compositing data, for the 
same parameter, from various data sets for the purpose 
of eliminating noise. 

3.47 Statics — Correction applied to seismic data to 
nullify the effect of elevation differences encountered 
along profiles, as well as the effect of a low velocity 
weathered layer. 

3.48 Suppressed Layer — Layer lacking a response 
because of its small thickness and/or contrast in 
physical property with the surrounding environment. 

3.49 Terrain Correction — Correction applied to 
measured gravity data to nullify the effect of irregular 
topographic relief in the immediate vicinity of the 
station of measurement. Charts are used to calculate 
the required correction. For local surveys in flat areas, 
this correction may not be required. 

3.50 Transition — Linear or exponential variation of 
a physical property with depth. 



3.51 Transverse Resistance — Product of the 
thickness and resistivity of a geoelectrical layer. 
Conventionally written as T= hxp (ohm.m 2 ). 

3.52 Two-Electrode (Pole-Pole) Configuration — 

Configuration in which one current and one potential 
electrode is kept at infinity (more than 10 times the 
distance between active electrodes) and perpendicular 
to the profile along which the other two active 
electrodes are moved. The geometric factor is 2na, 
where a is the distance between the active electrodes. 

3.53 Vibroseis — Seismic survey in which a vibrator 
is used as a non-destructive source instead of an 
explosive to generate controlled frequency seismic 
waves in the ground. 

3.54 Wenner Configuration — Collinear four- 
electrode configuration of potential and current 
electrodes in which all the electrodes are equidistant, 
that is, the separation between potential electrodes (a) 
is 1/3 rd the separation between current electrodes. The 
geometric factor is 2na. 

4 UNITS OF MEASUREMENT 

Table 1 lists the parameters and units of measurement 
in common use. 

5 PURPOSE OF SURFACE GEOPHYSICAL 
SURVEYS 

5.1 Surface geophysical surveys play a vital role in 
groundwater exploration. Surveys can be used to 
conduct either shallow subsurface investigation that 
may be needed for many environmental related projects 
or deeper investigations that may be required to identify 
productive aquifers. Also, surveys can be used to 
estimate the thickness of weathered zones, delineate 



Table 1 Commonly Used Geophysical Techniques and Units of Measurement 

(Clause 4) 



SI No. 


Method 


Technique 


Physical Property Involved 


Unit for Parameters Measured 


(1) 


(2) 


(3) 


(4) 


(5) 


i) 


Electrical 
resistivity 


Sounding 
Profiling 


Resistivity 


Ohm-m 


ii) 


Magnetic 




Mag. Susceptibility 
Mag. Field intensity 


Gammas 
Nano Tesla 


iii) 


Electromagnetic 


VLF 

HLEM 

TEM 


Conductivity /Resistivity 


Inphase/Quadrature Component ( percent ) 

-do- 

Secondary/Primary Magnetic Field (percent ) 

Voltage decay, Ohm-m, V Sec. 


iv) 


Seismic 


Refraction 


Wave velocity 


m/s 






Reflection 


Acoustic Impedance 


m/s 






(High Resolution) 






v) 


Induced polarization 




Chargeability 


Milli-second 


vi) 


Self potential 
(Electrokinetic) 




Natural Potential 


mV 


vii) 


Mise-a-la-masse 
(Charged body) 


Charged body 


Development of potential 


raV 


viii) 


Gravity 




Density (Lateral variation) 


Milli-galon 
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bed rock topography, demarcate fracture geometry, 
identify the presence of limestone cavities and/or 
paleochannels, and to assess quality of groundwater. 
Furthermore, surveys can be used to assess 
groundwater pollution and the movement of plumes, 
define vadose zone characteristics required for waste 
disposal or artificial recharge projects, demarcate sea 
water intrusion, differentiate between aquifers and 
aquitards, monitor the quality and direction of 
groundwater movement etc. Surface geophysical 
measurements are also used to estimate hydraulic 
parameters of aquifers. They are increasingly used 
because they are rapid and cost effective and they 
supplement direct methods such as drilling. 

Surface geophysical methods can be grouped into two 
categories — natural field methods and artificial source 
methods. Commonly used natural field methods 
include gravity, magnetic and self-potential methods 
which measure variations in earth's gravity field, 
magnetization and natural electric potential of rocks. 
Microgravity techniques, which detect changes in 
ground water storage, can be used to identify saturated 
cavernous limestone features. Artificial source methods 
measure the response of the subsurface to artificially 
induced energy like seismic and electromagnetic waves 
and electrical currents. These methods include 
electrical resistivity (see IS 15736), induced 
polarization, very low frequency (VLF) 
electromagnetic, controlled- source electromagnetic, 
seismic refraction (see IS 15681), and occasionally, 
seismic reflection. 

5.2 One of the well developed method is Ground 
Penetrating Radar (GPR) which is a high-resolution 
system for imaging subsurface using electromagnetic 
(EM) waves in the frequency band of lOHz-2000 MHz. 
It is used to detect the anomalous variations in the 
dielectric properties of the various subsurface materials. 

The GPR system consists of the following: 

a) A source for transmitting EM waves. 

b) Receiver for detecting EM waves reflected 
from different subsurface features. 

c) Control and display unit for synchronization 
between transmitter and receiver as well as 
recording, processing and display of data. 

5.2.1 Benefits of GPR 

a) Portability 

b) Application is non-destructive 

c) Rapid in data acquisition, and 

d) High-resolution subsurface imaging. 

5.2.2 Applications of GPR 

a) Detection of fracture zone, 



b) Determination of depth to water table, 

c) Location of sinkholes and cavities, 

d) Detection of anomalous seepage, and 

e) Mapping of archeological remnants. 

5.2.3 Limitations of GPR 

a) Penetration depth and ability to resolve targets 
at a depth is dependent upon the prevailing 
underground conditions. 

b) Highly conductive soils subsurface material 
render the GPR method ineffective. 

c) Sufficient electrical contrast between the 
target and the host materials is necessary. 

d) Interpretation of GPR data is subjective. 

6 PLANNING 

Surface geophysical surveys need to be carefully 
planned in order to meet project objectives. Planning 
should include the following aspects. 

6.1 General Considerations 

a) Effectiveness and accuracy of equipment and 
power supply, 

b) Easy operation and maintenance, 

c) Ready to use accessories, 

d) Suitability of vehicle for transportation, and 

e) Safety of equipment. 

6.2 Access to the Area 

a) Suitable access to the area/site, 

b) Permission to work in the area, 

c) Physical constraints in the area, 

d) Clearance along profile line(s), 

e) Noise and cultural disturbances, and 

f) Overhead power line. 

6.3 Equipment 

a) Maintenance should be performed as 
required; 

b) Should be stored in a stable, dust free, and 
dry environment; 

c) Pre-operation checking should be carried out; 

d) Power supply should be checked regularly; 

e) Precautions given for each equipment are to 
be observed; and 

f) Any deterioration in equipment condition 
should be rectified immediately. 

6.4 Safety and Precautions in Operation 

A safety code or plan should be developed prior to 
surveys to account for potential hazards in the field. 
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Common hazards include working with high voltage 
power lines, in electrical storms, in extremely remote 
areas, and with explosives. If possible, surveys should 
be conducted in dry weather periods to avoid damage 
to equipment by lightening. Unnecessary use of high 
voltage input should be avoided and care should be 
used when working with systems of 100 V or more, or 
with systems having 120 mA or more of current. In 
the event of rain or lightening, the current and potential 
cable connections should be removed from the 
instrument and no one should be allowed to touch the 
terminals. Even at a distance of 5 km to 6 km, 
lightening can damage the circuit. 

In seismic surveys, explosives should be handled by 
trained personnel stored safely. Overhead power lines 
should not be located near the shot hole, which should 
be dampened by water in the event that vibroseis or 
weight dropping is not used. Detonators should be 
always kept short circuited, even during transportation 
to the site. 

6.5 Planning of Survey 

Field crews should be informed of operational 
procedures prior to the survey. Profile lines should be 
straight and the distances between transmitter and 
receiver should be accurately determined. Spacings 
should be repeatedly checked or confirmed. Other 
considerations are itemized below: 

a) Crew should not touch the electrodes or the 
cable until instructed to do so by the operator. 

b) Movement of the crew near the profile should 
be restricted and the cable should not be 
passed through water or near high voltage 
power lines. Also, the crew should not stand 
in water in bare feet. 

c) Data should be plotted at the site so that errors 
can be removed or readings repeated. 

d) Electrodes should not be located near lateral 
inhomogeneity such as boulders in rocky 
terrain or buried objects such as pipe lines or 
telephone cables. 

e) Line should be checked regularly irrespective 
of the applied voltage. 

f) The charge (explosive) should not be placed 
in a highly-weathered zone so as not to overly 
dissipate the energy. 

g) For shallow investigations, the depth of 
weathering should be estimated by special 
shooting so that charge can be placed below 
the weathered zone. 

h) For EM equipment with multiple frequency 
selections, frequencies should be changed 
only after switching-off the instrument. 



j) In magnetic surveys ferrous objects should 
not be placed near the sensor. 

6.6 Quality Control in Field Data Collection 

Quality control considerations are a function of the 
selected equipment and the required level of accuracy. 
In any case, measurements should be repeated and 
profile orientations should be checked. 

6.7 Site/Area Details 

Investigators should become familiar with the local 
geologic and hydrogeologic characteristics of a 
targeted site prior to conducting a survey. 
Characteristics may include, but not be limited to, 
lineament details, lithostratigraphic information, water- 
level information, and water-quality information. A 
well inventory should be conducted to identify sources 
of pertinent data and information. 

Depending on the objectives of the survey, candidate 
sites for field surveys may be selected on the basis 
of existing information. Final site selection, however, 
should be based on a more rigorous study of 
geomorphic features and geological structures in the 
field. Local representatives may be consulted to help 
plan the surveys. Final site selection should be based 
on geophysical anomaly positions, accessibility, local 
conditions, and avoiding physical constraints such 
as electrical lines, metallic structures, crossing of 
roads, streams, or bridges, and topographic 
depressions. 

7 ELECTRICAL RESISTIVITY 

7.1 Purpose 

To identify groundwater-yielding zones (whether 
granular or fractured), zone geometry, variations in the 
chemical quality of groundwater, and the directions of 
groundwater movement (see IS 15736). 

7.2 Principles of Measurement 

A known amount of electrical current is first sent into 
the ground through a pair of electrodes. The potentials 
developed within the ground due to this current are 
then measured across another pair of electrodes on the 
ground. The distribution of current and equipotential 
lines in an electrically homogeneous subsurface is 
shown in Fig. 1 . The potential difference, V, between 
any pair of electrodes at the ground surface, P l P 2 , as 
shown in Fig. 2, is then calculated as 



V- 



2n\a b c d 



where p is the electrical resistivity of the 
homogeneous ground, / is the electric current with 
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Fig. 2 Different Electrode Configuration in Use in Electrical Resistivity Surveys 
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which the ground is energized, and a, b, c and d are 
the inter-electrode distances. Usually, both the current 
and potential pair of electrodes are placed in a straight 
line with the potential pair being placed inside the 
current pair to maintain a symmetry with respect to 
the inter-electrode distances. Two electrode arrays 
being used today are the Wenner and Schlumberger 
arrays as shown in Fig. 2. In the Wenner array, the 
electrodes are equally spaced while, in Schlumberger 
array, the potential electrodes are relatively close to 
one another as compared to the current electrodes. 
For the Wenner configuration of electrodes, the above 
equation becomes. 



2 a — 
I 



KR 



For Schlumberger configuration, it becomes 

2 



-MN 



MN 



V 

1 — KR 
I 



As shown in the above equations, when the resistance 
R is multiplied by K (a constant called the 'spacing or 
geometrical factor' which depends upon the spacing 
between current and potential electrodes), it gives the 
value of p, the resistivity of the ground. If the ground 
is homogeneous, the value of p gives the true resistivity 
of the medium or the ground. However, since the 
earth's subsurface is multilayered, the value of pa 
provides what is called the apparent resistivity value. 
Along with the electrode spacing, the apparent 
resistivity value is a function of the thicknesses and 
true resistivities of the individual layers, and deducing 
the true resistivity value of any individual layer is a 
difficult proposition. In practice, as the separation of 
the current electrodes is step-wise increased, the 



current penetrates and becomes more focused deeper 
into the ground. A plot between the current electrode 
separation and the resultant electrical resistivity value 
yields a curve known as vertical electrical sounding 
curve (in short VES). 

There are two ways of interpreting the VES data. The 
first involves matching the field curve with master 
curves that have been prepared for multi-layered 
system with different combinations of resistivity and 
thickness. The second method is computer aided where 
the VES curve is calculated for an initial 'best guess' 
model of the system and then adjusted by successive 
iterations to match observed curves. The matched 
model curve is assumed representative of a subsurface 
with the same layering and resistivity as indicated in 
master curve. 

In resistivity profiling, an electrode array (Wenner array 
is generally preferred) is moved in a line from one point 
to another to record variations in resistivity along a 
profile. The technique is helpful in locating lateral 
inhomogeneities owing to the presence of resistive or 
conductive bodies such as dykes, saline water bodies, 
etc. Significant resistivity contrasts occur between dry 
and water- saturated formations, and formations with 
fresh and brackish or saline water. Sands of various 
grain size, clays, weathered and fractured granites and 
gneisses, sandstones, cavernous limestones, vesicular 
basalts, etc, all have defined but overlapping ranges of 
resistivity. The resistivity ranges shown below for 
different materials are generalized and may vary 
significantly based on local hydrogeological 
conditions. 

7.3 Instruments 

A resistivity survey is carried out using an instrument 
known as a resistivity meter. These meters typically 
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employ either a direct current or a very low frequency 
alternating current type of suitable wattage and may 
also be equipped with noise filters and digital displays 
of current input and measured voltage. Measurement 
accuracies for many resistivity meters typically fall 
within the micro-volt range. Meters with multi- 
selection constant voltage or constant current input 
are desired. Required accessories include rugged 
winches/reels of insulated base with 200 m to 500 m 
of PVC insulated single conductor cable, multi-strand 
thin wires of low electrical resistance, a rechargeable 
or non-rechargeable direct-current power source, 
small diameter stainless steel rods/stakes and 
hammers, non-polarizing electrodes and connectors, 
hand-held walkie-talkie sets, and surveying 
equipment. 

7.4 Field Procedures 

There are a variety of electrode configurations used in 
resistivity surveys. The co-linear, symmetrical 
quadripole spread of the Schlumberger configuration 
for sounding and the Wenner configuration for profiling 
are the most popular. In the Schlumberger 
configuration, the practice is to move current electrodes 
outward while keeping the closely-spaced potential 
electrodes fixed at the center so long as a measurable 
potential difference is obtained. When the potential 
difference becomes so small that it cannot be accurately 
measured, the potential electrodes are expanded, 
always with the proviso that their separation does not 
exceed one-fifth that of the current electrode 
separation. Conventional sounding commences when 
the potential electrode spacing is equal to one-fifth of 
the current electrode spacing. Successive spacing of 
electrodes is usually increased in geometric 
progression, with each current electrode spacing being 
1.414 times the preceding one. As such there should 
be equal distribution of 6 points to 8 points in each log 
cycle of double log graph paper used for plotting the 
apparent resistivity curve. Spacing can be increased 
by 2 m to 5 m to study minor changes. In Wenner 
sounding, the potential electrode spacing is set at one- 
third the current electrode spacing through-out the 
survey (that is, all four electrodes are equidistant and 
moved outward for successive measurements). In 
hardrock areas, radial soundings (soundings taken at a 
site along 4 to 8 different directions) may be useful for 
studying fracture orientation and for correcting depth 
estimates. 

The Schlumberger and Wenner configurations each 
have advantages. The Schlumberger configuration 
requires less manpower and cable and, because 
electrode movement is relatively small, the effects of 
near surface lateral inhomogenities on the signal is 
minimized. Also, shifting of the curve with potential 



electrode changes are smoothed. The Wenner 
configuration has the advantage of giving higher 
potential values because the potential electrodes are 
equally spaced with the current electrodes. 

For sounding curves, apparent resistivity values are 
plotted against half current electrode separation for 
the Schlumberger configuration, against inter- 
electrode spacing for the Wenner configuration, and 
against the distance between the current and potential 
dipoles for the dipole-dipole configuration. For radial 
soundings polar diagrams are also prepared. In 
profiling with Wenner/Schlumberger/dipole-dipole 
electrode arrangements, the configuration (of fixed 
electrode distance) is moved along a straight-line 
profile taking measurements at fixed spacings (station 
intervals). In gradient profiling, current electrodes are 
planted well apart, say 800 m to 1 200 m, and the 
central one-third space is scanned by a potential dipole 
of 10 m to 20 m in length, at a station spacing of 5 m 
to 10 m. Gradient measurements can also be made 
along closely spaced (50 m apart) parallel profiles 
within the central one-third space without changing 
the positions of the more distant current electrodes. 
Groundwater flow and velocity can be measured using 
a rectangle configuration of potential electrodes 
placed midway between the two current electrodes in 
such a way that a uniform electric field exists near 
the potential electrodes. 

In profiling, apparent resistivity values are plotted 
against stations on arithmetic graph paper. The center 
of the potential electrode spacing is the point of 
measurement for the Wenner and gradient 
configurations. For the dipole-dipole configuration, 
the point of measurement is between the current and 
potential dipoles. When attempting to trace a fracture 
zone, because low resistive readings in a single profile 
may be erroneous and misleading, profiling should 
be taken along 2 to 3 parallel profiles located 50 m to 
100 m apart. Also, profiling should be preceded by 
test soundings to select optimum electrode spacings. 
At least one profile should be conducted with a small 
electrode spacing (5 m to 10 m) to understand the 
effects of near-surface resistivity variations on deeper 
information and to reduce ambiguities. In the Wenner 
configuration, the effects of near surface 
inhomogeneity can be reduced by an off-set 
arrangement of electrodes and by taking averages. 

Selecting a site for a survey should serve its purpose. 
In the event a geophysical anomaly is identified at a 
point which is not accessible for drilling, its extension 
should be identified by observing some parallel 
profiles. If site is near a concrete structure like road, 
building or bridge, profile should be laid in such a 
way that potential electrodes do not fall within 10 m 
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of the structure and are on homogeneous ground. Rock 
debris and building materials lying in vicinity of the 
electrodes should be removed. 

Locations of the current electrode positions should be 
identified before starting the survey. Electrode 
locations are accurately measured from the center and 
small pits/holes are made compatible to the size of 
potential electrode base and moisture condition of soil. 
In dry soil conditions, sufficient water should be put 
in pits/holes before placing the electrodes in the 
ground. Care should be taken to ensure that the 
electrodes are in proper contact with the ground. In 
case electrode location falls on dry, compact soil and 
sand, sufficient water is put in the hole by removing 
the electrode and placing again to minimize contact 
resistance at the electrode. For small electrode 
spacings, the electrodes should not be driven more 
than 40 mm to 50 mm to maintain it as a point 
electrode. For large spacings, the entire length of the 
current electrodes can be driven into the ground. 

The potential values should be higher than 5 mV and 
in no case below 1 mV. Minor variation in potential 
brings in noise in the apparent resistivity curve. Small 
potential values are generally obtained with large 
electrode spacings for which the "geometric factor" 
is quite large and relatively small inaccuracies with 
the geometric factor gives relatively large variations 
in apparent resistivity. 

Ideally, current circuit should not offer path resistance 
other than signal resistance. Therefore, in practice it 
should be ensured that cable resistance as well as 
contact resistance is minimum at both of the current 
electrodes. Contact-resistance can be reduced by 
driving the current electrodes deeper, and by putting 
saline water in electrode pits. If necessary, an 
additional one or two electrodes could be planted near 
the current electrode, about a meter apart, and 
connected in parallel to the main electrode. 
Alternatively, a sheet of tin foil placed in a watered 
pit can be a very effective current electrode. 

With the Schlumberger array, when potential electrode 
positions are changed, repeat measurements should 
be made for at least two of the earlier current electrode 
positions (with new potential electrode position) for 
overlapping curve segments. 

It is necessary to plot data during operation, so that 
trend of the curve is known and data points with noise 
can be repeated and also, the spacing to terminate 
measurements could be properly chosen (for instance, 
when bedrock is indicated by a steeply ascending 
portion of the VES curve). Accuracy of the data 
depends on the sensitivity of the instrument to measure 
potential differences, to filter out noise by stacking 



and displaying the standard deviation of measured 
values, and to correct measurement of electrode 
distances and their alignments. 

7.5 Processing of Data 

Sounding curves obtained by the Schulmberger 
configuration are generally discontinuous with upward 
or downward shifting of curve segments because of 
the shifting of potential electrodes. Shifting should be 
in a prescribed manner if there is no lateral 
inhomogeneity. Sounding curves can be smoothed by 
shifting the curve-segments up or down, depending on 
the type of curve (whether ascending or descending). 
Conventional shifting of the curve depends on the 
relative resistivities of the layer sequence. When the 
potential electrode spacing is increased, the depth of 
investigation is somewhat reduced, producing a curve 
that ascends upward and not downward. Difficulties 
involving the shifting of curve segments can be 
overcome by observing the trends of nearby soundings. 
Shifting of curve-segments could also be due to surface 
inhomogeneities near the potential electrodes. 

Surface inhomogeneities near the current electrodes can 
also be recognized by distortion in the sounding curve. 
A sharp curvature of the maximum value in the sounding 
curve is not indicative of a resistive layer of regional 
extent, but rather a lateral surface inhomogeneity. Curves 
that suddenly rise or fall with changes in the position of 
the current electrode indicate the presence of a 
lithological contact. In such areas, other nearby sounding 
curves can help smooth the distorted curve and identify 
which current electrode has caused the shifting. 

7.6 Interpretation 

Qualitative interpretation of sounding curves can be 
made visually to identify the type of curve and to 
demarcate areas with similar types of curves (for 
example ascending/descending type or H, A, K, or Q 
type curves for various combinations of multi-layered 
subsurface resistivity variations (Fig. 3). Quantitative 
interpretation of resistivity sounding data is based on 
empirical or semi-empirical methods in which the field 
curves are smoothed and matched with a variety of 2 
layer and 3 layer theoretical master curves along with 
the auxiliary point charts. This graphical technique 
involving a sequence of partial curve-matching where 
two or more homogeneous and isotropic (assumed) 
layers are combined in a single anisotropic (introduced) 
layer, which is equivalent to another fictitious single 
homogeneous and isotropic layer. Interpreting results 
from soundings made in relatively layers is difficult 
and to some extent depends on the skill and experience 
of the interpreter, and on the availability of local 
geological information. 
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Development of computer-based inversion techniques 
has greatly aided investigators in interpreting results. 
With these techniques, parameters values for targeted 
layers can be obtained from the iterative adjustment 
of estimated 'guessed' values to match the curves 
observed in the field. The equivalence and error 
analysis are done and also some of the layer 
parameters can be fixed (through borehole 
information), while inverting. A number of VES from 
an area can be interpreted simultaneously as 'batch 
interpretation' required for a regional consistency in 
results. In some of the inversion programmes, a guess 
model is not required. One such programme, which 
gives results with deeper layers in increasing order of 
thickness, is not very useful. In another, the curve is 
inverted by the process of "peeling-off" of layers. 
Here, the resistivity of the last layer is not correctly 
estimated because the process involves extrapolation 
of the last segment of the curve. Alternatively, it is 
advisable to interpret curves by forward modeling as 
well as automatic inversion as the former gives scope 
of incorporating geological information, while the 
latter provides more highly-resolved results, as well 
as an estimate of the error in final parameter values. 

An empirical approach is used in interpreting sounding 
curves from hard rock areas. For some of the current 
electrode spacings, input current becomes 
automatically high (when a constant current source is 
not used) and the curve shows a descending kink for 
that spacing. Statistical analysis has shown that a linear 
correlation exists between kinks observed in the curve 
and the depths of saturated fractures encountered in 
the borehole. The distance of the current electrode 
position for which a kink is observed is almost same 
as the depth to the fracture. 

Resistivity profiling data are interpreted qualitatively. 
From gradient profiling data, the ratio of the resistivity 
'low' (indicating saturated fracture zone) to the 
background 'high' is computed and calibrated with the 
borehole results, if available. That is, similar ratios in 
the same hydrogeological environment should indicate 
similar fracture zones. Besides the ratio of a low 
(anomaly) to background value, actual values as well 
as the steepness of the anomaly are also considered 
for an indication of the depth to the anomaly source. 
Quantitative interpretation should also include essential 
aspect of standardization of parameters through 
available borehole information. The interpretation is 
modified with the inflow of drilling data. 

7.7 Advantages 

The electrical resistivity method is cost-effective and 
employs non-destructive field techniques. It is effective 
in assessing the quality of ground water and therefore 
can be used to locate the saline/fresh ground water 



interface, or saline water pockets. Resistivity contrasts 
associated with presence or absence of ground water 
can be used to delineate the geometry of aquifers and 
zones favourable for ground water accumulation. This 
method also provides useful information on lithologic 
characterization, depth to resistive bedrock, direction 
of ground water flow, orientation of fracture zones, 
and the locations of faults and paleo-channels, as well 
as cavities in limestone. The method also can be used 
for specific environmental applications such as 
delineating the area and extent of ground water 
pollution, identifying zones suitable for artificial 
ground water recharge, soil salinity mapping, and 
reclamation of coastal saline aquifers (see IS 15736). 

7.8 Disadvantages 

Overlapping resistivity ranges and a very wide range 
of resistivity makes it difficult to characterize ground 
water targets by their resistivities unless standardized 
locally. Also, the accuracy and resolution of the 
response decreases with increasing depth and 
decreasing contrasts in resistivity. Finally, like other 
methods based on potential theory, is limited in its 
predictive application (see IS 15736). 

7.9 Limitations 

The presence of very high or very low resistivity surface 
soils can affect interpretation. While the former 
increases the contact resistance, the latter masks the 
signals coming from deeper layers. These presence of 
such soils can be problematic because they can 
attenuate a considerable percentage of the input signal 
going into the subsurface, as well as the output signal 
coming back from deeper zones. The resistivity low 
that may result from the presence of a conductive top 
soil/overburden may be mistaken for a suitable target. 
It is therefore essential that a profile with a very small 
electrode spacing is also conducted to identify the top 
soil conductivity effect. Cable resistance and contact 
resistance affect the ground resistance (measured 
signal) which is generally too low. 

Because the response of a resistivity profile is 
dependant on two parameters, that is, on the geometry 
and resistivity of the targeted layer, there is no unique 
solution and a number of equivalent models are found. 
While conducting soundings on a multi-layered earth, 
it is observed that the parameters of intermediate layers 
could be altered to a certain extent, keeping either the 
ratio of thickness-to-resistivity or the product of 
thickness and resistivity constant. This would not 
produce any appreciable/detectable change (within the 
accuracy of the observation) in the shape of the 
resistivity sounding curves. This phenomenon is known 
as equivalence, the effect of which is pronounced if 
the layers are thin. It cannot be resolved by a single 
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technique but requires the support of independent 
information for fixing either of the interpreted 
parameters or by obtaining the same parameters 
through joint interpretation with other techniques. The 
response is dependent on the depth and resistivity 
contrast of the target. Thin layers or layers with less 
resistivity contrast with the surrounding are 
suppressed. 

In a layered sequence, the interfaces between 
successive layers having monotonously increasing or 
decreasing orders of resistivity cannot be distinguished 
accurately, particularly at greater depths, because of 
the transitions in resistivity. 

Dipping layers distort the measurements and produce 
ambiguities. The presence of inhomogeneties, either 
at the potential or current electrodes, produces distorted 
or shifted curves which can be difficult to interpret. 
Induction of pseudo-anisotropy due to repetition of 
resistive and conductive layers results in an error in 
depth estimation, requiring corrections by calculating 
the coefficient of anisotropy through transverse 
resistance and longitudinal conductance. 

8 SELF POTENTIAL 

8.1 Purpose 

To obtain information on the direction of natural 
ground water movement or seepage, or the movement 
of ground water induced by pumping of a water well. 

8.2 Principles of Measurement 

Measuring natural self potential is one of the oldest 
and simplest methods in geophysics. The natural 
electrical potential [also known as Self Potential or 
Spontaneous Potential (SP)] has two components, 
namely the electro-kinetic component and the 
electrochemical component. The electro-kinetic 
component or streaming potential component of SP is 
linked to the flow of ground water, making its use 
effective in ground water exploration. Streaming 
potentials increase in the direction of groundwater flow 
and the gradient of anomaly is related to the magnitude 
of flow. That is, a map of equal SP values would reflect 
the direction and magnitude of flow. The magnitude 
of streaming potential is generally low, being of the 
order of millivolts. 

8.3 Instrument 

Because the amplitudes of the anomalies produced by 
the streaming potential can be quite low, a 
potentiometer or resistivity meter capable of measuring 
in the millivolt range is required. Also, a micro- 
processor based stacking facility would help rejection 
of noise. 



8.4 Field Procedures 

An SP survey is carried out along 20 m to 50 m spaced 
parallel profile lines or along radial lines originating 
from a borehole in 8 to 12 directions. Station intervals 
can be kept at 2 m to 10 m, depending on the objective 
of the survey. Water-filled electrode pits are 
constructed in advance of the survey so that potentials 
are stabilized. If possible, inhomogeneities located near 
the potential electrode should be removed while 
making the pits, or the pits should be constructed an 
adequate distance away from inhomogeneities. While 
taking measurements, the presence of inhomogeneities 
are to be recorded. The electrodes should be firmly 
placed into the pit. If porous pot-type potential 
electrodes are used, they are usually kept connected in 
a tub containing a copper sulphate solution for 8 h to 
12 h prior to their use in order to minimize potential 
differences due to the electrodes themselves. 

There are two techniques used in SP field 
measurements — the total field measurement 
technique and the gradient or leap-frog measurement 
technique. In the total field measurement technique, 
one of the potential electrodes is kept fixed as a base 
or "reference electrode" at a site geologically suitable 
(that is without much variation in potential), while the 
other electrode is moved along the profile lines. If the 
reference electrode is shifted, a new reference electrode 
is tied in with the previous one and measurements are 
overlapped. In the gradient measurement technique, 
both electrodes are moved along profiles lines with a 
fixed separation. The distance between the electrodes 
is kept very small. The total field measurement 
technique is preferred, as it gives large values of 
potential difference and the error associated with 
electrode polarization is less in comparison to that in 
the gradient configuration. All measurements should 
be completed in the minimum time possible to avoid 
drift due to polarization. 

To study the direction of ground water movement 
induced by pumping, measurements are usually taken 
around the well before pumping and then repeated after 
pumping for a reasonable duration, after switching- 
off the pump. 

8.5 Processing of Data 

In the total field measurement technique, the data are 
reduced to a common point and corrected for drift. 
Polarization effects are reduced by linearly 
interpolating it between the measurements. Corrected 
data can be plotted as profiles or as iso-potential 
contour maps. 

8.6 Interpretation 

Interpretation of SP data can be difficult because a 
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number of factors, including noise, may distort the 
streaming potential anomaly. Also, the order of 
magnitude of the noise may be same as that of the 
anomaly (± 5 mV to 10 mV). Qualitative interpretation 
is done by studying the amplitude and wavelength of 
the anomaly and then matching the anomaly with 
patterns for known source geometries. Smooth 
anomalies with longer wavelengths indicate a deeper 
source mechanism. Shorter wavelengths and higher 
amplitudes indicate a shallower source. Vertical 
lithologic contacts or structural discontinuities give 
steep, asymmetrical anomalies with an amplitude 
dependent on the resistivity ratio. Streaming potentials 
are reduced with increasing clay content. Interpretation 
of SP profile data is often facilitated by comparison 
with geoelectrical and/or geological sections and 
topographic profiles. 

8.7 Advantages 

Self potential is a relatively inexpensive method for 
obtaining useful information on the lateral as well as 
the vertical movement of ground water flow or 
seepage. In areas of polluted groundwater, flow 
induced by pumping can be traced to plan preventive 
measures. The method is also useful for locating 
shallow water-filled cavities in limestones with 
appreciable ground water flow. 

8.8 Disadvantages 

Interpretation of results may be difficult because SP 
anomalies are often complex and of low amplitude. 

8.9 Limitations 

In long profile lines, SP anomalies are affected by 
telluric current variations which may be of much higher 
order. Also, SP anomalies are affected by magnetic 
storms and may yield erroneous anomalies on sloping 
ground. The presence of near surface inhomogeneities, 
conductive overburden, variations in soil moisture, 
electrochemical effects, conductive/ resistive bodies in 
the subsurface, overhead power lines, and corroded 
pipe lines will obscure the anomalies due to the 
streaming potential. Measurements are also affected 
by the location of the reference electrode and the 
watering of electrodes during measurement. 

9 FREQUENCY DOMAIN ELECTROMAGNETIC 
(HORIZONTAL LOOP) 

9.1 Purpose 

To delineate saturated fracture zones in hard rocks and 
to estimate the thickness of weathered zones. 

9.2 Principles of Measurement 

In conventional electromagnetic (EM) surveys a 



transmitter radiates electromagnetic waves (primary 
field) that penetrate the ground. When the primary field 
encounters a conductor, that is, a body of limited 
extensions with an electrical conductivity higher than 
its surroundings, eddy currents are produced in the 
conductor. A secondary electromagnetic field (in a 
direction opposite to the primary at the conductor) is 
produced by the eddy currents and the resultant field 
is measured by a receiver, placed at a given distance, 
in the form of in-phase and quadrature components 
[see Fig. 4A]. The receiver also measures the primary 
field. The resultant field is either measured as a 
percentage of the primary field or its direction relative 
to the vertical is recorded. The magnitude, direction 
and phase angle (which is the time delay of the resultant 
field in relation to the primary field) of the resultant 
field can be used to locate a conductive body and obtain 
its parameters. There are several ways to conduct EM 
surveys by varying the position and orientation of 
receiver and transmitter loops, namely, vertical loop, 
horizontal loop, Turam, etc [see Fig. 4B]. Overall, in 
EM exploration it is generally assumed that there exists 
a conductivity variation in the subsurface and that the 
conductive target is located within a non-conducting 
(resistive) surrounding, or that the conductivity of the 
target is much higher than the surrounds. 

The EM method has advantages over the resistivity 
method in that the latter has difficulties in sending a 
current through a highly resistive surface layer, such 
as those often found in deserts or in compact rocky 
terrains. Also, the change in penetration depth can be 
obtained by changing the frequency of the transmitted 
electromagnetic wave, as well as the transmitter- 
receiver coil separation. Because anomalies on 
groundwater targets in hard rocks are caused by 
conductivity contrasts between the saturated zone and 
the surrounding dry medium, a higher contrast would 
provide a better response. The electromagnetic method 
has been used widely in groundwater exploration, 
occasionally to compliment the resistivity method and 
help resolve ambiguities in interpretation. 

A commonly used technique is the Horizontal Loop 
Electromagnetic (HLEM) method, also known as the 
Slingram method. HLEM surveys are controlled-source 
surveys in which the transmitter can be operated at a 
number of frequencies and transmitter-receiver coil 
separations. The transmitter and receiver coils are 
placed in the same horizontal plane. HLEM profiling 
with a number of frequencies and transmitter-receiver 
separations gives a depth-wise distribution of electrical 
conductivity. That is, a reduction in the frequency of 
EM waves and/or an increase in the transmitter-receiver 
separation would provide deeper information. As in 
resistivity surveys, a conductive overburden of varying 
thickness can create a problem in quantitative 
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interpretation and in detecting the target. The primary 
(incoming) field suffers attenuation etc and the depth 
to the target is sometimes overestimated as in the case 
of the electrical resistivity method. 

In-phase and quadrature components of resultant 
magnetic field expressed as percentage of primary 
magnetic field are measured [see Fig. 5A, 5B]. 
Resultant field is a function of conductivity, frequency 
and coil separation. Hence, measured values depend 
on the response parameter a = U. a CO L 2 , where u. is 
magnetic permeability, a is ground conductivity, CO 
is angular frequency and L is transmitter receiver coil 



separation. In electromagnetic surveys the term 
'conductivity' is preferred as response is generally 
proportional to conductivity. 

9.3 Instrument 

The instrument is comprised of a transmitter, a receiver, 
and the console. The transmitter can be operated at a 
number of high and low frequencies, usually in the 
range of 100 Hz to 10 000 Hz. The instrument should 
have repeatability of readings. Teflon-coated cables 
of 50 m, 100 m and 200 m lengths are often used for 
connecting the transmitter to the receiver. 
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9.4 Field Procedures 

HLEM surveys are usually conducted as profiling in 
combination with electrical resistivity surveys. Even 
a few resistivity soundings could be conducted in a 
study area to define the geoelectrical layering prior 
to detailed HLEM profiling. 

Profile lines are laid across the probable strike of the 
target conductor. Profile lengths are kept much longer 
than the expected lateral extent (geometry) of the 
anomaly. 

Station intervals are kept 10 m to 25 m apart, 
depending on the objective of the investigation and 
likely target dimensions. Spacings between the 
transmitter and the receiver should be accurately 
measured. For a spacing of 100 m, a maximum error 
of 200 mm is permissible. 

Transmitter and receiver coils should be placed on 
the ground horizontally and properly oriented towards 
the receiver (that is, in the same plane), corrected 
through inclinometer or tiltmeter. The receiver and 
transmitter should be properly aligned. To measure 
the response, the transmitter and receiver coils are 
moved in unison to successive stations, keeping inter- 
coil spacing fixed. In-phase and quadrature 
components of the resultant field can then be measured 
at available frequencies. Changes in frequency are 
indicative of various depths of penetration. Therefore 
each station with in-phase and quadrature data 
measured at 4 frequencies to 8 frequencies represents 
an EM depth sounding. Measured values represent 
the information obtained from the center of the 
transmitter-receiver coil separation. Profiling may be 
conducted at two or more separations of receiver- 
transmitter coils. 

9.5 Processing of Data 

Data can be plotted for inphase and quadrature 
components together or separately for different 
frequencies. Noise in the data can be eliminated by 
visual inspection. Phasor diagrams can be prepared 
to estimate the layer parameters. 

9.6 Interpretations 

Interpretation of the target anomaly can be done 
qualitatively as well as quantitatively. The width of 
the anomaly is equal to sum of the thickness of the 
conductor and the coil separation. Quantitative 
interpretation includes curve matching with available 
or generated theoretical models that have been 
previously developed for various subsurface 
conductivity distributions, depth-to-thickness ratios, 
and conductor altitudes. 



The presence of conductive overburden increases the 
amplitude of the anomaly. At higher frequencies, the 
quadrature component response produces a base level 
shift and may reverse or become negative. The 
conductor appears to be buried deeper and more 
conductive. 

Sounding data can be presented as a phasor diagram 
and interpreted with available sets of such diagrams 
that have been prepared for various layered earth 
models. The presence of a conductive surface layers 
rotates the phasor diagram clockwise. 

Using an initial guess model and certain assumptions, 
the sounding data can be inverted by software to give 
layer models at each point. Interpretation becomes 
more useful if some borehole information is available 
to identify the character of geologic structures 
producing the response. 

9.7 Advantages 

EM field operations are fast and cost-effective and 
can produce voluminous data. The instrument can be 
operated at a number of frequencies and coil 
separations for depthwise information. There is no 
need of ground (galvanic) contact, so no operating 
problem of current injection or of contact resistance 
in areas of highly resistive surface layer and also no 
noise introduced in the data because of near surface 
inhomogeneities. The EM method provides better 
lateral resolution and assessment of rock conductivity 
than does the electrical resistivity method. 

The method requires less coil separation for deeper 
information than do resistivity soundings. As a rule 
of thumb, penetration depths for HLEM are 1.5 times 
the transmitter-receiver coil separation distance, 
compared with a maximum penetration of about one 
quarter of the current electrode separation required 
of the Schlumberger resistivity sounding. 
Consequently, given a favourable subsurface 
conductivity distribution, much deeper information 
can be obtained by covering less ground space. Also, 
multi-frequency data give deeper information, that is, 
depth of penetration is not constrained by coil/ 
electrode spacing as in the resistivity method. 

9.8 Disadvantages 

Success of the method depends on getting a 
meaningful interpretation of the data, which in turn 
depends on the conductivity characteristics of the 
overburden through which primary field penetrates 
and returns, introducing two phase lags. That is, the 
HLEM technique is preferred to detect a conductive 
target through less conductive overburden, which may 
not be always available. 
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For layered earth interpretations, models are highly 
simplified which may not be the true condition. 
Detection of deeper layer is difficult. Skin-depth plays 
a significant role, as depth of exploration depends on 
relative conductivity of deeper layers. 

In phase component is affected by topographic 
variations and it is always essential to standardize the 
response of the model through known borehole 
information in the area. 

9.9 Limitations 

The problem of equivalence exists. Presence of 
conductive overburden or surface layer induces phase 
lag and ambiguity. It is difficult to differentiate 
anomalies due to overburden variation and those due 
to variation within the bedrock. 

Interpretation of layer model is not very accurate for 
highly conductive and resistive surface layer, that is, 
for the high contrast in conductivities. 

10 TRANSIENT (TIME DOMAIN) 
ELECTROMAGNETIC 

10.1 Purpose 

To delineate aquifer zones in a conductive surrounding 
and delineate conductive saline ground water zones. 

10.2 Principles of Measurement 

The transient electromagnetic (TEM) method is 
related to the frequency domain (continuous wave) 
electromagnetic methods by the Fourier transform. 
Instead of making measurement at different 
frequencies, in TEM methods the decay of an induced 
voltage is measured at a number of sampling times. 
A constant current is passed through the transmitter 
loop which produces a static primary magnetic field. 
When the transmitter current is abruptly switched off 
the static magnetic field decays and, due to the 
associated flux changes, currents are induced in 
conductors in the ground. This current flowing in a 
horizontal closed path below the transmitter loop 
produces a secondary magnetic field. The change in 
amplitude of secondary magnetic field with time 
induces a voltage in the receiver coil. Response 
normalized by the primary field is measured at 
selected time intervals after switching-off the primary 
field. Because response depends on resistivity of the 
ground, measurements can yield geoelectrical 
characteristics of the ground. Immediately after 
switching off, that is at early time stage induced 
current is concentrated near the surface of the earth. 
Since, the maximum amplitude of induced current 
diffuses downward and outward, deeper geoelectrical 
information can be obtained as time increases, that 
is, at later stages. The transient field decays quite fast. 



The shape of transient curve (voltage decay versus 
square root of time or apparent resistivity versus square 
root of time) does not represent depth- wise resistivity 
variations as it could be assessed from conventional 
dc apparent resistivity curve. Actually, the depth of 
exploration is a function of time (and current flowing 
in the transmitter loop) and does not depend on 
transmitter-receiver separation. 

10.3 Instrument 

Transient electromagnetic system comprises a receiver 
and a transmitter loop unit. Transmitter loops of 
different sizes are used for exploring different depth 
ranges. TEM instrument uses constant current 
waveform consisting of equal periods of time-on and 
time-off. A variety of TEM equipments are available 
with stacking facility. The TEM measurements are 
made in a time range of 6 [is to Is after switching off 
the primary current. The latest measurement time is 
determined by level of noise. For shallow groundwater 
exploration measurement up to 10 millisecond to 
30 millisecond is done. 

10.4 Field Procedures 

The technique can be employed for sounding as well 
as profiling. For profiling moving transmitter- 
receiver configuration is used. Three types of 
transmitter-receiver configurations are employed in 
TEM soundings, namely, grounded line, central loop 
and loop-loop configurations. The grounded line 
configuration is used for deep soundings, while 
central loop and loop-loop configurations are used 
in shallower applications. The dimension of the 
transmitter loop in central loop configuration 
depends on the depth to be explored and is selected 
based on field testing. 

The transmitter loop dimensions range between about 
30 m x 30 m to 500 m x 500 m to explore shallow 
zones to depths of about 2 500 m. For better resolution 
at early time a small loop size is desired. Large loop 
size at later times provides better signal. A peak 
current of 2 A could be sent through the loop of 30 m 
x 30 m for shallow exploration. Higher amperage 
(20 A) and large loop size is used for deeper 
exploration. The receiver measurements can start at 
6 |is after switching-off and therefore shallow zones 
can be investigated. The latest time could be up to 10 
millisecond to 30 millisecond depending on the level 
of noise. The minimum detectable signal ranges from 
10" 6 V/A.m 2 to 10 12 V/A.m 2 . A group of four to six 
persons are required as crew. 

10.5 Processing of Data 

The voltage decay versus time observed data are 
converted to apparent resistivity versus time data. 
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10.6 Interpretation 

Induced voltage decay curve does not present direct 
picture of the subsurface geoelectrical condition as in 
case of electrical resistivity. Data are normalized for 
transmitter and receiver parameters and converted to 
apparent resistivity. Apparent resistivity versus square 
root of time is plotted on double log graph paper. 
Curves are interpreted either by curve matching or by 
software packages for inversion and forward modeling 
. Problem of equivalence exists in this technique also. 

10.7 Advantages 

Data scattering is not observed in central loop 
sounding as in case of electrical resistivity sounding. 
It is least affected by lateral variations in resistivity 
as the induced current flows in rings around the 
receiver and also transmitter loop size is not changed 
frequently. Resolution is high in shallow central loop 
soundings. It has better resolving capability for S {hip) 
equivalence and can be used with other techniques 
that respond better to resistive layers (that is direct- 
current electrical sounding) to help resolve ambiguity. 
Compared to electrical resistivity sounding, smaller 
area/smaller loop size is required for survey to achieve 
same order of depth. Thus, it can be conducted easily 
in confined areas. To probe deeper, transients at later 
times are recorded. It is highly sensitive to 
conductivity changes, that is, a highly conductive layer 
underlying conductive clay overburden is detected 
better than a resistive layer. 

10.8 Disadvantages 

TEM equipment is quite expensive. Practically, to 
overcome noise, transmitter loop size has to be 
increased to investigate deeper targets. 

A good estimate of first and the last layer may not be 
possible, due to equipment constraints. To get 
information for near surface layer very early stage time 
data is required. 

Target of limited lateral extents may not give a good 
match in inversion (due to 3-D effects). 

Resistive fresh water aquifers underlying thick clay 
overburden may not get detected. Also, if the first layer 
is quite thick and resistive, the deeper relatively 
conductive layer may no get detected. 

10.9 Limitations 

Thin resistive layers can not be detected. Transient EM 
noise at later time stage restricts the length of time 
during which transient can be sampled and thus deeper 
information cannot be obtained unless transmitter loop 
size or primary current flow is increased. Transient 
sounding for deeper exploration requires a large area 



for loop layout compared with the straight strip 
required for co-linear electrical resistivity arrays. 
However, for similar depths of exploration, dc sounding 
methods sample a much larger volume of ground and 
the data are therefore more susceptible to 
inhomogenities and reduced resolution. 

Information on ambient noise at the measurement 
location is necessary. 

Technique may not be useful if resistivity-thickness 
contrast is comparable with measurement uncertainty. 

11 VERY LOW FREQUENCY (VLF) ELECTRO- 
MAGNETIC 

11.1 Purpose 

To delineate conductive water bearing fracture zones 
in resistive hard rock and to determine approximate 
thickness of overburden. 

11.2 Principles of Measurement 

VLF method is a type of electromagnetic method in 
which only receiver is in control of the operator. 
Transmitters are fixed stations located at great distances 
(up to several thousand kilometre) from the survey area. 
There are several such transmitting stations around the 
world, which are continuously emitting 
electromagnetic waves in frequency range of 15 kHz 
to 30 kHz for navigation purposes. Though the term 
VLF indicates very low frequency, the technique uses 
quite a high frequency for geophysical applications. 

At large distances from the transmitter, radiated EM 
waves travel into the ground as plane wave with a 
horizontal magnetic and electric field and a vertical 
electric field all mutually perpendicular. These plane 
waves (primary field) penetrate the earth surface and 
in case a conductor (relatively conductive saturated 
fracture zones) is present eddy currents are created in 
it. A secondary field with arbitrary orientation is 
generated due to the current induced. The resultant 
magnetic and electric fields are not in phase with the 
primary and so are elliptically polarized. In VLF 
surveys, secondary field due to eddy currents is 
measured by a sensitive receiver. 

The technique has directional advantages as well as 
limitations. Saturated fracture zones in hard rocks 
oriented in-line with transmitter are picked up with 
relative accuracy. VLF has a better resolving power 
because of higher frequencies, but is effective in 
detecting shallow fracture zones only. The high 
frequency radiations are attenuated fast with depth. 
Also, in presence of conductive over-burden, 
attenuation is fast and technique becomes ineffective 
in detecting deeper fracture zones. Thus, VLF response 
becomes very susceptible to unwarranted near surface 
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inhomogeneities and apparently presents a plethora of 
anomalies. VLF receivers can measure tilt of the major 
axis of polarization ellipse and the ratio of minor to 
major axis (known as ellipticity). VLF receivers can 
measure in-phase part, which is approximately equal 
to tilt and the out-of-phase (quadrature or imaginary) 
part, which is approximately equal to the ellipticity of 
vertical component of secondary field expressed as a 
percentage of the horizontal primary field. Electric field 
normal to primary magnetic field is measured to 
measure apparent resistivities. 

11.3 Instrument 

Several types of instruments are available at present. 
Besides instrument, conventional surveying 
accessories are required to lay out profiles along 
desired orientations. 

11.4 Field Procedures 

VLF method is easy to operate in field. A transmitter 
of strong and clear transmission is selected. If there is 
option of selecting 2 or 3 transmitting stations, they 
are to be selected, keeping orientation of fracture zone 
in mind, to get the maximum response. 

Parallel profile lines are laid perpendicular to the 
direction of transmitting station that is, along the 
direction of primary field. 

Profiles are laid 25 m to 100 m apart and station interval 
is kept at 10 m. Length of profile is generally kept large, 
more than a kilometer, to study the regional trend. 
Selection of station interval becomes quite significant 
where anomalies show high rates of curvature. 

Orientation of receiver with respect to transmitter is 
adjusted by a method given for the instrument selected. 

For some instruments, there is no need of keeping 
receiver at specific orientation with respect to the 
transmitter. 

Apparent resistivity of surface layer can also be 
measured by some of the instruments, using two sensors 
connected to instrument and placed on ground about 
5 m to 10 m apart at each station along the profile. 

Operational procedure varies with type of instrument. 
In some of the instruments data are direct, digitally 
displayed, while in other they are recorded by obtaining 
a minimum intensity of sound signal adjusting the 
instrument in various positions/orientations. 

Position of transmitting station with respect to the 
movement of operator, that is, to his right or to his 
left is to be noted for interpretation of the "cross- 
overs" of inphase and quadrature components. 



Accuracy of data depends on signal-to-noise ratio and 
selection of transmitter with reference to the orientation 
of target. 

In field operation, repeatability of readings is to be 
ensured. Instruments in which minimum sound is 
observed, accuracy may vary with the operator and 
affect the readings. 

11.5 Processing of Data 

In-phase and quadrature components are plotted along 
profile line. Noise in the profile is identified and 
removed. If data on parallel profiles are available 
contour maps can be prepared for in-phase and 
quadrature components. 

11.6 Interpretation 

Data can be interpreted qualitatively as well as 
quantitatively. Being a reconnaissance survey method 
it is mostly used for qualitative assessment. Anomalies 
are identified and interpretation of depth and lateral 
extents of targets and their conductivities are assessed. 
Mostly, the technique is used to demarcate lateral extent 
of a target. 

Anomalies being affected by the presence of thick 
conductive overburden, assumption and simplifications 
are required in interpretation. Effect of conductive 
overburden and conductive host rock surrounding 
should be studied in detail from the available literature, 
before making any inference. 

Quantitative interpretation is also attempted, in which 
experience of interpreter plays a significant role. In 
highly resistive terrain, ratio of in-phase to out-of-phase 
response is proportional to conductivity of the target. 

For quantitative interpretation data can be filtered using 
Fraser and Karous-Hjelt filters. 

Fraser filter is used for in-phase data, which show cross- 
over response. It turns cross-overs into 'peaks and 
troughs' and reduces sharp noise. Filter enhances those 
anomalies, which resemble its shape. 

Karous-Hjelt filter is used to determine the subsurface 
distribution of current, which is responsible for the 
measured magnetic field. Current density pseudo depth 
sections are obtained for the purpose. 

Quantitative interpretations can also be attempted for 
layered earth model. 

11.7 Advantages 

VLF survey is fast and economical in field operation and 
used for reconnaissance in delineating saturated fracture 
zones in hard rocks. Surveys can be made in areas where 
surface layer is highly resistive and high contact resistance 
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would be encountered in galvanic resistivity surveys. 
Lateral disposition of conductive zone is delineated 
accurately. It gives fast estimation of surface soil/ 
overburden resistivity. Use of higher frequency range, 
enhances resolving power in detecting closely spaced 
conductivity discontinuities. Detectability of target 
increases in resistive surrounding. 

11.8 Disadvantages 

The amplitude (or even the continual presence) of the 
VLF primary field cannot be always guaranteed at 
the receiver. 

To get proper response and detection, there is 
restriction on orientation of the target zones. 

Depth of investigation depends on resistivity of the 
surrounding media, and is drastically reduced if 
surface layer is highly conductive. 

Instrument is expensive and may not deliver as much 
information of the subsurface as resistivity method 
can, except for reconnaissance. 

Because of high frequency used, measurements are 
highly susceptible to variation in resistivity and 
thickness of overburden. Most of the anomalies are 
generated by the variations in overburden alone and 
can be mistaken for the underlying fracture zones. 
So, the data profiles are noisy. Data obtained is a 
function of operational procedure and hence 
ambiguities are induced. 

11.9 Limitations 

Because of high frequency, the fields are attenuated 
and phase shifted. 

Conductivity resolution is effective over a frequency 
range. 

Secondary field attenuates fast and skin-depth is small 
in highly conductive formations. 

In conductive terrain maximum depth of penetration 
is half skin-depth for the medium surrounding the 
target or overlying it. 

12 SEISMIC REFRACTION 

12.1 General 

Seismic refraction technique is quite useful to map 
areas suitable for existence of potential aquifers. It is 
used to determine thickness and differentiate 
compactness of sediments, subsurface layering, 
delineate weathered zone thickness, bedrock 
topography, identify fracture zones and palaeo- 
channels. Sometimes the technique is very effective 
in differentiating saturated and un-saturated zones 
(see IS 15681). 



12.2 Principles of Measurement 

Seismic vibrations are created artificially on the surface 
of the earth either by explosion of dynamite (high 
energy), impact at the ground surface either by heavy 
and accelerated weight drop ( medium energy) or by 
hammer (low energy). Vibrations thus created spread 
to underground spherically in all directions and their 
arrival at different distances at the surface of the earth 
are detected by sensors, planted on the ground, known 
as Geophones. The responses of the geophones are 
recorded in seismograph with timer circuit so that times 
of arrival of these vibrations, called seismic waves, 
from the shot point, where the vibrations are created 
to the detector points, where geophones are planted at 
different distances, are accurately measured ( in milli 
second). The greater the compactness of the medium, 
the higher the velocity. 

In sedimentary or loose alluvial formation, the velocity 
of seismic wave propagation increases if the medium 
gets saturated with water. Similarly seismic wave 
velocity in weathered rock will be conspicuously less 
compared to compact rock system. Of the various 
seismic waves, generated, Longitudinal wave, also 
known as Primary wave or, in short P wave, is fastest 
and first to be detected. Thus, in refraction seismic 
work conducted for ground water exploration, 
propagation of only the P wave through different 
subsurface layers is considered. 

The subsurface consists of different layers and is not 
homogeneous. The compactness of the layers 
generally increases with depth and as a result, the 
deeper layers are expected to have seismic wave 
velocity greater than that in the overlying material. 
This condition, which is necessary for the refraction 
method to be successfully applied, creates refraction 
of the down moving seismic wave, incident to the 
interface of the two layers at a particular angle, called 
critical angle. At the interface, the refracted wave, 
sometimes called head wave, moves with velocity of 
the lower layer. As a result, after some distance 
between the source and the receiving geophones, the 
refracted wave takes over the direct wave and is first 
to reach the detector. 

To explain Snell's law here, highly relevant to 
refraction principle, if incident ray enters the first 
medium with P wave velocity Vj at an angle a with 
vertical and emerges as refracted wave in the second 
medium with P wave velocity as V 2 at an angle (3 with 
the vertical [Fig. 6A], it may be proved based on simple 
principle of optics that 

sin a V, 



sin P V 2 
When angle a and velocity contrast between two media, 
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that is, difference between V l and V 2 becomes such that 
P becomes 90°, the above equation simplifies to 



sin a 



sin 90° V, 



or 



sin a ■■ 



V, 



(Necessary condition for critical refraction is V X <V 2 ) 

This is the condition of total refraction and refracted 
wave grazes the interface and every point of the 
interface becomes active and according to Huygens 
Principle, acts as secondary source to generate waves 
that is transmitted back to the first medium as head 
wave. Angle a is called critical angle. 

In seismic refraction technique, head wave 
propagation and detection at the ground surface are 
of main concern. The principle as explained may be 
seen in Fig. 6B. The wave front analysis is the best 
way to understand the refraction phenomenon and an 
illustration is shown below to appreciate Huygen's 
principle for the generation of head wave. Schematic 
presentation of refraction seismic survey and plot of 
travel-time curve is shown in Fig. 6C. An impression 
of 24 channel seismograph is presented in Fig. 6D 
where actual impression of direct and head wave 
arrivals at detectors (geophones), planted at different 
distances from the shot can be seen. 



In the schematic Fig. 6(A) of travel time curve, shown 
above, f ic is the intercept time and X ci is the critical 
distance, where direct and refracted wave reaches at 
the same time. 

Plot of time of arrival of P wave versus shot-detector 
distance makes it possible to determine P wave 
velocity in different layers (medium) and either from 
the value of f ic or X ci depth to interface or thickness of 
first layer can be determined with following formulae 



H-- 



X : V.-V 



2 vv 2 +v, 



The calculation of depth from intercept time will be 

as 



H 



t 



1 _ 1 



Value of velocity indicates the compactness, degree 
of water saturation in layers, degree of weathering 
and fracturing in rocks, etc. The depth to different 
subsurface layers and seismic velocities in them 
(multi-layered system) can be calculated in the same 
way as indicated above for two layered earth. There 
are simple mathematical derivations with which dip 
of different layers also can be calculated. Reciprocal 
time method techniques have been adopted by various 




6A Snell's Law 
Fig. 6 Technique Measurement by Seismic Refraction — Continued 
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Fig. 6 Technique Measurement by Seismic Refraction — Continued 
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Fig. 6 Technique Meosurement by Seismic Refraction 



authors based on time distance discrete analysis 
techniques to find out depth of interface at different 
detector points and these are highly relevant to 
accurately decipher discontinuity in the refractor 
interface, like presence of fractures, faults, etc. 
Approximate plus-minus technique in refraction 



survey, based on wave front techniques are also 
espe cially effective in mapping interface 
discontinuities, especially in solving hidden layer 
problems. 

Finally, a tentative range of seismic P wave velocities 
in various lithologies is shown below as: 



SI 


Formation 


PWave 


No. 




Velocity 
(m/second) 


(1) 


(2) 


(3) 


i) 


Loose alluvial sediments 
(dry/partially saturated) 


300-800 


ii) 


Loose and saturated alluvial 
sediments 


1 200-1 500 


iii) 


Sandstone 


1 500-3 000 


iv) 


Quartzitic sandsone 


3 000-5 000 


v) 


Limestone 


3 000-6 000 


vi) 


Cavernous limestone 


2 000-3 500 


vii) 


Deeply weathered rock 
(Granite/Charnockite/ Khondalite) 


1 000-2 000 


viii) 


Partially fractured granite 


2 000-3 500 


ix) 


Compact granite/gneiss 


4 500-6 500 


x) 


Charnockites 


3 000-4 500 


xi) 


Khondalites 


2 000-3 000 



Velocity Range m/second 



100-500 500-1 000 

(4) (5) 



1 000-2 500 2500-5 000 >5000 
(6) (7) (8) 
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12.3 Instruments 

Seismic data acquisition system comprises geophones, 
amplification, filtering and recording system. For 
seismic refraction surveys, geophones of natural 
frequency of about 10 Hz are used. The amplification, 
filtering and recording constitute the instrument, 
which could be analog or digital. It can have a 
continuous record of all the seismic signal reaching 
the geophones or only record a strong signal. 

12 geophone to 24 geophone channels are employed 
to study several hundred metres of depth and for very 
shallow depths one geophone channel is sufficient. 

The seismic sources are hammers, weight drops and 
explosive. The hammer and weight drops are used for 
shallow exploration. 

12.4 Field Procedure 

Surveys are carried out along straight profiles laid by 
surveying. In refraction surveys total length of profile 
will be much greater (generally 3 to 5 times) than the 
depth to be investigated. A series of 12 geophones or 
24 geophones are planted at regular interval along 
the profile. 

The equal spacing up to 10 m between the geophones 
is kept for shallow investigation and 30 m to 50 m for 
investigating depth of 300 m to 400 m. Surveys are 
to be carried out in the areas of minimum cultural 
noise and vibrations, or during the period when the 
noise becomes minimum, for example during nights 
in the urban areas. Surveys are to be avoided along 
roads with traffic, railway lines and during high winds, 
etc. 

Profiles are shot from both the ends to ascertain the 
dip of the interfaces and for reciprocal time 
information. The shot holes are made by a small 
drilling system or weight dropping is applied to 
generate the energy. 

12.5 Processing of Data 

Differences in elevation of shot point and geophone 
is corrected by bringing them to a common datum by 
subtracting or adding the travel time, if the velocity 
is known. The presence of near surface low velocity 
layers of varying thickness will produce an apparent 
variation in the depth to the underlying interface. 

12.6 Interpretation 

A preliminary interpretation is done in the field after 
a shot is fired. There are a number of quantitative 
interpretation methods, out of which the intercept time 
method and plus-minus methods are commonly 
applied. The intercept time method is based on Snell's 



law. Travel time curves available for the shots from 
both ends of the profile are analyzed to check the 
reliability and dip of the layers. The plus-minus method 
of Hagedoorn is an approximate method but fast and 
simple, based on the reconstruction of wave front. The 
intercept-time and plus-minus methods can be 
combined for practical solutions. The detailed 
literatures on all the quantitative interpretation 
methods are available. Advanced computer softwares 
like Seisrefa or Q Seis, are available to directly process 
digitized data obtained from seismographs. Very 
sophisticated seismographs with digital record 
facilities are available where various kinds of signal 
processing techniques are available which earlier was 
not possible with analog recorders. 

12.7 Advantages 

A stratigraphic sequence could be established, if the 
stratigraphic breaks are associated with velocity 
contrast. Accurate depth can be estimated as a support 
to electrical resistivity interpretations. It has better 
resolution in vertical section. Bedrock topography and 
fracture zones can be deciphered with reasonable 
accuracy. 

12.8 Disadvantages 

Field operations are not easy. Cannot be operated in 
noisy area or urban areas. Use of explosives prevents 
its general application. Also, the overall cost of surveys 
(equipment, field operation, field arrangements and 
explosive etc) sometimes prohibits its application for 
ground water exploration. 

12.9 Limitations 

The assumptions made are highly simplified and 
seldom met with. The effect of lateral seismic velocity 
changes within the layer will be an overestimation or 
underestimation of depth. The variation in seismic 
velocity in horizontal direction from vertical direction 
will cause an overestimation of thickness, if not 
considered. 

Presence of 'low velocity layer' — the blind zone — 
contradicts the necessary assumption for seismic 
refraction. It can cause overestimation of depth to the 
layers underlying the low velocity layer. Presence of 
'thin layer' (thin with respect to its depth of 
occurrence) induces error in depth estimation. 
However, both of these conditions may be detected 
and compensated for in some modern methods of 
interpretation such as Palmer's Generalized 
Reciprocal Method. The presence of dipping layers 
and their effect will induce error in estimation of 
velocity and therefore profiles are shot from both the 
ends for correction. 
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13 SEISMIC REFLECTION 

13.1 Purpose 

To delineate subsurface layering, channels and thin 
layers, geological structures and bedrock topography, 
etc. 

13.2 Principles of Measurement 

Seismic reflection in its simplest form uses a near 
surface source to generate seismic wave which travels 
through the subsurface formations and get reflected at 
interfaces (reflectors) having a density-velocity 
contrast. The velocity of seismic wave and the time 
taken to reach the reflector and return to the surface 
yield the depth to the reflector after putting some 
corrections to the data. In seismic reflection, for an 
angle of incidence less than the critical angle, the 
seismic ray is reflected back towards the ground surface 
and is received by a large number of receivers or 
detectors known as geophones. 

13.3 Instrument 

A variety of seismic reflection equipment are available 
for measurements in standard frequency range. For 
high frequency-high resolution seismics special types 
of equipment and receivers are used. Several types of 
energy sources are available, namely, explosive, 
vibroseis, weight dropping, etc. For explosive source 
a shot hole is required to place it below the water table. 
Best quality of data is obtained by it. Vibroseis 
generates vibrations in an oscillatory manner for a small 
duration instead of impulsive one. Weight dropping is 
a low energy and less expensive source. Non-explosive 
sources are useful in populated area. 

Receivers also known as geophones convert ground 
motion into electrical signal. The signal from each 
geophone station is amplified and passed through cable 
to the instrument. For high resolution surveys, high 
frequency geophones are used. 

13.4 Field Procedures 

13.4.1 Sources 

A number of seismic source are used for shallow 
seismic reflection work. Hence the choice of seismic 
sources for shallow investigations is made on the 
considerations of the cost, repeatability, convenience, 
efficiency, safety, soil and ground conditions and safety. 
The sources which are used are as given below: 

a) Sledge hammer, 

b) Weight drop, 

c) Blasting caps and explosives, 

d) Rifles and shot guns, and 

e) Wackers (Minisosie). 



The blasting cap has a significant frequency range of 
100 to 200 Hz. The rifle bullet impact (depending 
upon calibre, bullet velocity and depth of penetration) 
has a frequency range of 100 to 200 Hz. Rifles and 
Wackers (Minisosie) have frequency range of 80 to 
100 Hz and explosive hammers and weight drop of 
50 to 200 Hz. 

13.4.2 Receivers (Geophones) 

Geophones having natural frequency of 50 to 100 Hz 
are commonly recommended for shallow seismic 
reflection work. One factor that effects geophone — 
ground coupling. It is therefore recommended that 
the size of spike or base plate of the geophones should 
be increased. The geophones are planted in the ground 
at regular intervals and are connected to the 
seismograph. The layout of the geophone is known 
as geophone array. A geophone array may have one 
or many geophones having different lengths and 
geometry. Array greater than 1 to 1.5 m is designed 
to attenuate the surface wave noise which adversely 
affect the reflection signal at frequency above 150 
Hz and array of 1 to 1.5 m are too small to attenuate 
effectively the surface wave lengths in excess of 10 
m. Therefore, the use of array in shallow seismic work 
is rather limited and single geophones are commonly 
used. 

13.4.3 Recorders (Seismographs) 

Until recently, engineering seismograph used only 
8 or 10 bit. Analog/Digital (A/D) conversion, they 
had a dynamic range not exceeding 60 db for a 10 bit 
A/D conversion. There was a need for addition of 
sufficient dynamic range, and use of filters especially 
low cut filter for carrying out seismic reflection work. 
The use of filters and signal enhanced seismic is 
essentially to increase the signal and noise ratio. 
Unlike seismic reflection methods normally applied 
in deep hydrocarbon exploration, its use in shallow 
projects require use of high frequency. It is therefore 
essential that seismograph use a high sampling 
interval, preferably 0.4 msec to 0.2 msec. Some of 
the new generation seismograph are 24-48 channels 
and uses 16-24 bit A/D. 

13.4.4 Seismic Survey Set Up 

The following two procedures are in vogue in shallow 
seismic reflection: 

a) Common depth point (CDP), and 

b) Common optimum offset technique (COF). 

In seismic reflection survey it is possible to record a 
number of shots such that a reflection point is 
common. This setup is called a common reflection 
point (CRP), a common depth point (CDP), or a 
common mid point (CMP). The degree of multiplicity 
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of data from a particular location is known as CDP 
fold. Most present day recordings use atleast 6 folds. 

13.4.5 Field Parameters 

The field parameters to be considered in designing a 
seismic reflection survey programme include purpose, 
a target depth, record wave length, sampling intervals, 
seismic source, geophone characteristics, source 
geophone distance, spread length, gain and filter setting 
on the seismograph. To obtain all these, it is therefore 
necessary that test seismic survey to be first carried out 
prior to actual survey to decide on optimizing the 
parameters mentioned above. Comparison of deep and 
shallow seismic reflection parameters are given in 
Table 2. 

Table 2 Equipment and Field Specifications for a 
Typical Survey 



SI 


Equipment 


Field Specifications 


No. 






i) 


Seismograph 


24 or more channels 

1 00 db minimum dynamic range 

Flat frequency response from 50 to 

1 000 Hz 

Wide choice of analogue filter settings 

Large spatial digital data storage 

(>10Mb) 


ii) 


Geophone 


50-100 Hz, 7.5 cm spike 


in) 


Source 


Sledge hammer or Weight drop or 
Buffalo gun, 12 gauge 


iv) 


Computer 


Pentium IV processor or better 
80 GB hard disc 
12 V power supply 




Record length 


0.5 s 




Sample interval 


0.4 or 0.2 ms 



Record length 


0.5 s 


Sample interval 


0.4 or 0.2 ms 


Pre A/D filter 




v) Geometry 




Spread type 


On end spread 


No. of channels 


24 minimum 


Geophone spacing 


2 m or 1 m 


Geophone/station 


1 


Shot spacing 


2 m or 1 m 


Record length 


0.500 s 


Nominal fold 


12 


Nominal CMP 


1 m 


spacing 




vi) Data Processing 




Gains and polarity 


Gain removal, Spreading gain corr 


Field statics 


Variable velocity 




Single layer straight 


Velocity analysis 


Spectra every 50 CMPs 


Residual statics 


Single pass, pilot 




From selected traces 


Stack 


Mean stack, selected 


Filter 


Bandpass 




Low cut 30 Hz, High cut 300 Hz, 




Notch 60 Hz as applicable 


NOTE — There can be many variations of the same survey. 



13.5 Acquisition and Processing of Data 

Acquisition of seismic reflection data requires 
experience and skill. It involves a number of steps 
including amplification and filtering. Data processing 
is done to improve the quality of data, that is, to 



enhance the signal to noise ratio. Software is available 
to process the data in a sequence. Processing 
parameters include resampling, filtering, 
determination of field statics, stacking, deconvolution 
and migration, etc. The end product is a seismic time- 
section. 

13.6 Interpretation 

Interpretation involves velocity determination and 
time to depth conversion and identification of 
reflectors. Useful interpretation of seismic-time 
section includes identification of litho-facies, 
structures and inferences on hydrogeological 
conditions. Deducing hydrogeological inferences and 
stratigraphic boundaries from seismo-geological depth 
section needs integration with geological information 
and depositional environment of the sediments. 
Computer softwares is available for interpretation of 
data. 

13.7 Advantages 

It is most accurate method to determine depth to the 
lithological interfaces and locate geological structures 
accurately. Thin layers can be detected by high 
resolution seismics. 

13.8 Disadvantages 

Field operation and processing of data are highly 
expensive. A large number of crew is required for the 
surveys. It needs a great deal of planning and local/ 
legal clearance. In using source as weight dropping a 
high level of noise is generated. Inappropriate selection 
of processing parameter may give wrong information. 
Though seismic reflection techniques have better 
accuracy and resolving power, because of high cost of 
seismic equipment and field operation, they do not 
enjoy wide applications in groundwater exploration. 

13.9 Limitations 

It cannot usually be operated in areas where noise, 
cultural disturbance or wind velocity is high. During 
field operation weather should be good and wind 
velocity should be low. In high-resolution seismic 
reflection surveys, that is high frequency surveys to 
identify thin layers the weathered layer (overburden) 
of varying thickness attenuates the amplitude, and for 
this source as well as the receivers are to be buried. 
Complexities in velocity distribution may indicate 
structures incorrectly. 

14 MAGNETIC 

14.1 Purpose 

To delineate subsurface structural features associated 
with contrasting magnetic susceptibilities, namely 



29 



IS 15897 : 2011 



basic dikes, quartz reefs, shear zones etc, in hard rocks 
that may contain useful supplies of groundwater or 
influence its flow. Also, it is employed for regional 
reconnaissance to narrow down the zone of interest, 
reducing cost of labour intensive detailed geophysical 
surveys. Magnetic surveys may also reveal shallow 
paleo-channels that can form useful aquifers. 

14.2 Principles of Measurement 

Variation in magnetic susceptibility of rocks causes 
variation in magnetic field intensity. Measurement is 
made of local variation in earth's magnetic field 
intensity, direction and gradient the magnetic 
anomalies caused by lateral variation in magnetic 
susceptibility of subsurface formation having 
magnetic mineral of different amount and nature. 
Measurement of magnetic field intensity along profiles 
is primarily a reconnaissance method (see Fig. 7). Unit 
of magnetic field intensity used in magnetic survey is 
nanotesla (nT) or Gamma. Intensity of earth's 
magnetic field ranges from about 25 000 nT at the 
magnetic equator to 65 000 nT at poles. 

Rocks containing varying quantity of ferromagnetic 
minerals may acquire both induced and permanent 
(remanent) magnetization. Intensity of magnetization 
is the resultant of these two components. 

Magnetic anomaly depends on the magnetic properties 
and geometry of the causative body, and the intensity 
and direction of earth's magnetic field at the point of 
measurement. Amplitude of anomaly ranges from as 
low as a few nT to thousands of nT. The geometry 
(wavelengths) of anomalies ranges from a few meters 
to hundreds of meters. In ground water exploration 
mostly local anomalies are significant. 

14.3 Instrument 

Proton precession magnetometer is mostly used to 
measure total magnetic field at 1 nT accuracy. It is a 
lightweight, portable equipment having a sensor and 
a console (see Fig. 8). It gives measured values in 
5-digit or 6-digit digital display. Instrument should 
have repeatability of readings. It should have multifold 
non-ferrous stands for holding the sensor at different 
heights. A total field intensity map for approximation 
of regional value should be available. 

14.4 Field Procedures 

Magnetic method is an easy field-operation 
geophysical survey for groundwater exploration. It is 
conducted mostly for reconnaissance and to support 
other methods, involving minimum time and money. 

Values of total magnetic field intensity are not 
measured at isolated spots but they are measured along 
selected profile line(s) at regular station interval for a 



meaningful interpretation of variations in the field 
intensity. 

Spacing of profiles, their orientations and station 
interval for measurements are decided based on 
objective, details required, size and depth of the body 
or anticipated size and amplitude of the anomaly and 
availability of space. Closer the spacing, greater the 
chance of detecting anomalies and accuracy in 
interpretation. Spacing and interval should be such that 
anomaly is not lost. 

Reconnaissance survey may be at large profile spacing 
and station interval while detailing should be done at 
close-grid. In ground water surveys profile spacing 
could be 100 m to 200 m and station interval 10 m to 
25 m or 50 m. Profile lengths should be 500 m to 
1 km. 

Profiles should be run normal to geologic strike which 
direction placed is known from local geological 
information. 

The base station has to be selected in an area of low 
magnetic gradient and where it is not likely to be 
disturbed for repeat base station measurements that are 
required for correction. 

Before starting the survey, a few readings should be 
taken at a particular station, to determine if there is 
any temporal magnetic disturbance (that is magnetic 
storm). If so, the survey should be delayed until the 
disturbance has subsided. 

Profile should not pass within about 100 m of buildings, 
power or telephone lines, buried pipe lines, near 
standing vehicles, railway tracks, fencing (metallic) 
which are extraneous source of magnetism. 

Measurements are to be repeated at selected base 
station by returning to it at convenient intervals (2 h to 
3 h) throughout the survey, depending on the length of 
profile and duration lapsed. 

For large-scale surveys where repeat measurements at 
one base station become cumbersome, more than one 
base station may be established. 

Time of observing the readings should be recorded 
frequently. Diurnal variation curve (plot of base 
readings with time) should be constructed for each 
day and corrections to be linearly distributed among 
the readings. Diurnal correction will be more effective 
if the variation is recorded almost continuously by 
another instrument deployed at the base station, as this 
correction is for external field and nothing to do with 
the instrument, deployed. 

Measurements are to be made by keeping sensor at a 
height of 2 m to 3 m above ground. Measurements by 
keeping sensor on the ground will be full of noise 



30 



IS 15897 : 2011 




MAGN E TIC \A LUE ( GAMMA) 

+ 




GROUND SURFACE 



BODY MAGNETIZED*INDUCTION IN THE EARTH'S FIELD 
(Q) EXAMPLE OF SKETCHED TOTAL FIELD ANOMALY 



GROUND SURFACE 





' H0R\-Z-0H7fi?f\ 



N 



BODY MAGNETIZED BY INDUCTION IN THE EARTH'S FIELD 
(b) EXAMPLE OF SKETCHED VERTICAL FIELD ANOMALY 



Fig. 7 Magnetic Anomalies 



31 



IS 15897 : 2011 




Fig. 8 Proton Precession Magnetometer 



because of the presence of varied material lying near 
it. 

Ideally some measurements of local rocks should be 
made with a susceptibility meter to assist with data 
interpretation. 

Accuracy of data depends on following factors: 

Instruments with 1 nT sensitivity will give sufficiently 
accurate data for ground water exploration purpose. 
Repeatability of reading is essential. Orientation of 
instrument has a role on the value measured. Reliable 
anomaly of interest is to be obtained by separating it 
from 'near' surface features. 

14.5 Processing of Data 

Magnetic measurements are corrected to eliminate 
variations and effects, which are not due to the body 
and have external source, and identify anomalies 
which may be characteristic of the subsurface 
formation or body. Correction or normalization is quite 
simple. 

Corrections made are mainly for : (a) diurnal, and (b) 
geomagnetic variations. Diurnal (solar and lunar) 
variations are caused by ionospheric influences and 
of periodicity about a day. Diurnal variations are of 
two types, namely the smooth, low amplitude (25 nT) 
regular variation and 'disturbed day' rapid, high 
amplitude (1 000 nT) irregular variation caused by 
magnetic storms, during which measurements are 



discontinued. 

'Quite day' diurnal variation is corrected and adjusted 
by recording variation through the day at one or two 
fixed base stations in the survey area and correcting 
all other measured values by interpolating/distributing 
the variation linearly over them. Closing or loop error 
in successive base readings of the order of 5 nT to 
10 nT is also distributed amongst the readings. Diurnal 
correction should be made for a chosen reference time. 

Geomagnetic variation relates the magnetic field 
variations with geographic coordinates of the 
observation point. For surveys covering a small area 
this correction is not required, but information is 
essential for interpretation. 

14.6 Interpretation 

Observed (corrected) values are plotted at stations 
along profile lines on a map and a contour map is 
prepared, showing lines of equal magnetic field 
intensity or different profiles of corrected values may 
be drawn. 

Data are interpreted qualitatively in terms of geologic 
structures and their lateral extents by visual inspection 
of contour trend; gradient of convergence or 
divergence and shape of the anomaly under 
consideration. 

A tentative geological interpretation for the area is 
made. 
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Quantitative interpretation is quite complex because 
of inherent ambiguity, equivalence and assumption 
made for simplification of model. 

In case of total field measurements anomaly 
component changes with latitude, for example, 
horizontal component at equator and vertical 
component at pole. 

Softwares is available for interpretations to iteratively 
match the observed anomaly with calculated anomaly 
for a defined simple geometric form and attitude. 

Closer the body to ground surface, steeper would be 
the gradient of magnetic anomaly and would be of 
smaller extents. Uncertainties in interpretation 
increase with complex geometry. 

14.7 Advantages 

Field survey is economical and coverage is fast. In 
hard rocks it is an essential reconnaissance technique 
to delineate concealed, shallow, geological feature 
favorable or unfavorable to ground water occurrences. 
Approximation of simple geometrical shape and size 
is possible. 

14.8 Disadvantages 

Equivalence problem (inverse potential) in 
interpretation exists. Anomaly is a complex response 
of the body involving several parameters. A number 
of assumptions are made for interpretation. It is difficult 
to operate in urban areas or areas with reinforced 
concrete structures. No inference is obtained on ground 
water quality conditions. 

14.9 Limitations 

Sensor cannot be operated near power lines. It may 
not be possible to continue measurements during mid- 
day. Measurements are to be discontinued during 
magnetic storms. Sensor cannot be kept on ground 
surface for measurements or near ferrous object. 

15 GRAVITY 

15.1 Purpose 

To ascertain change in ground water storage in valley 
fill aquifers, to locate dry or water filled cavities in 
limestone, to map bedrock topography and subsurface 
structures, where sufficient lateral density contrast 
exists. 

15.2 Principles of Measurement 

Gravity method involves measurement of minor 
variations in gravitational field of the earth caused by 
lateral changes in density of near surface rocks. Like 
magnetics, gravity method is also a natural source 



method (passive). It can be used as a reconnaissance 
method of survey dealing mainly with the relative 
variations. Gravity anomaly is expressed in milligal 
(gal being derived from Galileo), 1 gal = lcm/s 2 or 
1 milligal = 10~ 3 gal. Mean value of gravity at the earth's 
surface is 9.8 m 2 /s or 980 gals. 

Density of rock does not vary over a wide range. It 
varies from about 1.5 g/cm 3 to 3.5 g/cm 3 . Dry alluvium 
has density around 1 .5 g/cm 3 while water saturated 
alluvium has a density around 2 g/cm 3 . Density of 
sandstone and limestone varies from 2 g/cm 3 to 
2.8 g/cm 3 . While density of igneous rock varies from 
2.5 g/cm 3 to 3.3 g/cm 3 . Densities of different rocks 
overlap. 

15.3 Instrument 

To measure the relative variations a gravity meter is 
used. Basically, it is an extremely sensitive balance. 
Available instrument is lightweight, and transportable. 
Instruments with accuracy of 1 microgal are available 
for micro gravity measurements. 

15.4 Field Procedure 

The manner in which gravity survey stations are laid 
out and their spacing, depend on the type of 
investigation. It could be conducted along a series of 
profile lines across expected strike direction of 
subsurface feature. Measurement can be made at small 
station interval (5 m to 10 m) for detailed study to locate 
cavities in limestone terrain. Otherwise it can be 
conducted at a regular (selected) station interval of 
100 m to 500 m. Additional measurement stations are 
located in areas where changes in gradient are fast. In 
general, it is quite fast and economical in field 
operation. 

The measurement stations are preferably kept in flat 
area without any marked irregularity. Gravity base 
stations are conveniently located. Measurements are 
closed at the base stations. If required, more than one 
base station can be established. Topographic elevation 
of the stations are to be precisely measured and this is 
the most time-taking and expensive component of field 
operations. Placing of gravimeter is made by leveling 
to check its verticality. Time of measurements must be 
noted. 

15.5 Processing of Data 

A number of corrections are made to deduce the gravity 
anomaly due to rock density contrast alone. They are 
instrumental drift, latitude, elevation and terrain 
corrections. Elevation correction comprises free-air 
correction and Bouguer correction. Bouguer anomaly 
is obtained after applying these corrections to the 
measured gravity values. 
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15.6 Micro-gravity Measurements 

With high refinements achieved in measuring variation 
in earth's gravitational field to microgal range, it is 
possible now to resolve problems like determining 
change in ground water storage in a closed valley fill, 
identify subsurface cavities in limestone which may 
act as a potential ground water conduit in Karstic 
terrain. The method is becoming widely used in 
engineering investigations to detect natural and man- 
made cavities. Time-varying differential micro-gravity 
is a powerful technique for detecting changes in mass 
distribution beneath an area, like change in ground 
water storage. 

On earth the acceleration due to gravity varies between 
9.78 ms 2 to 9.83 ms~ 2 . These units are much too large 
and unit gal is preferred to identify structural features 
and change in subsurface density (an alternative unit 
in use today is the 'gravity unit' which is one-tenth of 
a milligal). Its unit is 1 gal = 0.01ms" 2 . Microgal is a 
further reduction of this unit and is defined as 
1 micro gal = 10" 8 ms 2 which are in use for search of 
voids, cavities and change in ground water storage 
which give anomalies of the range of few tens to 
hundreds of micro gals. 

The micro-gravity survey technique is a very speedy, 
accurate and cost effective geophysical technique for 
identifying natural cavities in limestone and other 
(man-made ) cavities, and monitoring changes in 
ground water storage in a closed well defined valley. 
The principle of the technique is to locate areas of 
contrasting density in the subsurface by collecting 
surface measurements of the lateral variation in the 
earth's gravitational field. 

Albeit simple to understand, a highly sophisticated data 
acquisition technique, maintaining the quality of data, 
advanced data interpretation are essential pre requisite 
in micro-gravity surveys. 

Whether the cavity is man-made or natural, it must 
have a way in and usually a way out otherwise there 
will no outlet by which material may be removed. 
Consequently, the 'effective' size of the target is 
dependent not only on the strict volume of the cavity, 
but also its connectivity, the secondary effects imposed 
by the cavity on the surrounding rocks which arise from 
the genesis of that cavity in its host rock. Several 
literatures are available where successful use of micro- 
gravity surveys has been made for engineering 
problems, detection of subsurface cavities and ascertain 
difference in ground water storage. Presently a wide 
scale application of micro gravity surveys are being 
made in ground water survey and exploration by USGS 
in western arid tract of USA in valley fills to ascertain 
change in ground water storage. 



15.7 Interpretation 

In gravity surveys, interest being in local anomalies 
(residuals) regional effect is eliminated from the 
Bouguer anomaly. Gravity effect of shallow and 
localized bodies could be enhanced and isolated from 
regional trend by second derivative analysis. 
Subsurface structures are inferred from maps of 
Bouguer anomaly contours, residual anomaly and 
second vertical derivatives. Gravity anomalies are 
positive in case of excess mass, that is, shallow bedrock, 
anticline, horst or ridges, etc, while negative anomalies 
indicate deficit of mass, that is, deeper bedrock, 
syncline, graben and cavity in limestones that is shape 
of the anomaly manifests the trend of subsurface 
geological structure. Gravity interpretation is 
ambiguous, as identical anomaly could be obtained 
for a variety of subsurface body placed at different 
depths (the inverse potential problem). When the 
causative body is shallow, the anomaly is sharp. 

Gravity anomaly is interpreted by indirect and direct 
methods. Anomaly could be matched and adjusted by 
iterations with the theoretical one for a model of known 
density contrast and geometry and depth of burial. 
Otherwise the anomaly is interpreted directly for shape 
and depth of the body. Accuracy and uniqueness in 
interpretation depends on reliable density estimates of 
the rock and information on geological and 
hydrogeological conditions. A number of computer 
software are available for interpretation. 

15.8 Advantages 

It is a fast reconnaissance method of survey and it is 
economical in field operation. Detailed micro-gravity 
mapping can precisely locate the cavities in limestone 
and buried channels. Detailed topography is required. 

15.9 Disadvantages 

Ambiguity in interpretation exists due to equivalence 
in interpreting potential field. Resolution is poor. 
Several corrections are required prior to the actual 
interpretation. 

15.10 Limitations 

Survey cannot be conducted in areas with cultural 
disturbances (that is vibrations). Precise elevation of 
measurement stations is essential. 

16 OTHER TECHNIQUES 

16.1 Induced Polarization 

16.1.1 Purpose 

To differentiate conductive clay mixed sand layers 
from conductive saline water-saturated zones. 
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16.1.2 Principles of Measurement 

On connecting two grounded current electrodes to a 
battery, a voltage between two points on ground can 
be detected. If current is switched-off, voltage does 
not become zero immediately but takes a small 
duration. The residual voltage also known as over- 
voltage is because of induced polarization effect. It is 
mainly concerned with surface polarization of metallic 
minerals induced by electric current and redistribution 
of positive and negative ions in the ground. Process 
of redistribution of ion can be classified into two 
groups: electrode polarization (exploited for metallic 
mineral exploration) and membrane polarization 
(useful for groundwater polarization). 

Membrane polarization occurs due to the presence of 
clay particles partially blocking ionic solution path. 
Surface of clay particle being negatively charged, 
attracts free positive ions from electrolyte and a double 
layer of charge is formed at the surface. Positively 
charged layer repels other positively charged ions and 
act as impeding membrane. On applying electrical 
potential, positive membrane is disrupted and positive 
charge carriers easily pass through while negative 
charge carrier accumulates. When current is switched 
off, redistribution of ions takes place and manifested 
as a decay of voltage between the two electrodes, that 
is, induced polarization effect is mainly the diffusion 
of ions. In either case, polarization is a surface effect 
and therefore greater the surface area of mineral particle 
or clay stronger would be the effect. In groundwater 
exploration membrane polarization is effective and silty 
formation shows maximum polarization. 

Measurements of induced polarization may be in either 
the time domain or frequency domain. The physical 
property and parameters measured are chargeability, 
frequency effect, metal conduction factor and apparent 
resistivity. The measurements are made for decaying 
voltage and the phase shift between receiver voltage 
and transmitter current. In time-domain, the measured 
quantity at the receiver electrode pair — the 
chargeability (m) — is the total area under the voltage 
decay curve. In frequency domain the change in 
resistivity at two different frequencies is measured. 
Frequency effect (f e ) is the difference in apparent 
resistivities at two frequencies normalized by higher 
frequency apparent resistivity. Metal conduction factor 
is calculated by normalizing the percent frequency 
effect by low frequency apparent resistivity. 

16.1.3 Instruments 

There are two types of instruments, time-domain and 
frequency domain. In time domain, a square pulse of 
current is passed through grounded current electrodes 
and the decaying voltage is measured between non- 



polarising potential dipoles at preset time intervals 
during current-off period. The output voltage is 
integrated over an interval and divided by the voltage 
applied to the current electrodes. Amplitude and 
duration of pulse are selected in such a way that decay 
of voltage is over, a second before the reverse pulse is 
transmitted. The standard cycle time is 8 with each 
on-off period as 2s. 

While in frequency domain, change in resistivity is 
assessed by measuring voltage at different frequencies. 
Modern IP receivers work in both frequency as well 
as time domain mode and also provide multi-dipole 
simultaneous reading facility. Frequency domain IP 
receivers are designed to measure both amplitude and 
phase over a broad band of frequency (2~ 14 Hz to 
2 +16 Hz). Instruments are microprocessor based 
with storage, facilitating fast recording of data, 
stacking of signals, filtering of noise and digital display. 

16.1.4 Field Procedures 

For groundwater exploration, IP sounding as well as 
profiling can be conducted. A variety of electrode 
configuration is used in IP survey. For routine field 
work Schlumberger, Wenner, two-electrode, dipole- 
dipole, pole-dipole and gradient configurations can be 
used. While pole-dipole and dipole-dipole 
configuration have low electromagnetic coupling, 
gradient configuration has high electromagnetic 
coupling. The survey with gradient array is relatively 
fast. 

In time domain surveys two parameters are selected. 
They are the time-period for which transmitter is 
switched on and off and the sampling interval of decay 
voltage. 

16.1.5 Processing of Data 

Normally, data are plotted at the mid point of potential 
electrodes or at the centre point of electrode 
configuration. When dipole-dipole electrode 
configuration with fixed inter-electrode distance is 
moved along a profile, data are plotted for each position 
of configuration at the intersection of two lines drawn 
towards the centre at 45°. This is repeated for each 
placement of configuration and for increasing 
separation between current and potential dipoles. 

16.1.6 Interpretation 

In saturated sediments IP effect is observed only when 
sands are mixed with some amount of clays. Coarser 
material show less effect than finer one. The effect is 
prominent in silts. Clean sand will not show any IP 
effect. A large quantity of clay segregation would show 
less polarization than the same amount of clay if 
dispersed. Polarization depends on the type of clay (that 
is polarization is more in montmorillonite than 
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kaolinite) but the variation is complex. Membrane 
polarization of clays decreases with increase in salinity 
of water. For fresh water saturated sediments, 
membrane polarization can also be correlated with 
hydraulic conductivity. 

Like apparent resistivity sounding curves, IP sounding 
curves can be interpreted by available theoretical 
sounding curve. Qualitatively it can be interpreted that 
pure clay will have low resistivity whereas sand will 
show higher resistivities. High resistivity associated 
with high chargeability would reflect presence of thick 
sand layers inter-layered with thin clays. For saline 
water zone chargeability may increase or decrease 
depending on the salinity. 

16.1.7 Advantages 

Induced polarization survey for groundwater 
exploration is taken up as supplementary to resistivity 
surveys. When it becomes difficult to identify the 
conductive clay mixed zones from the conductive saline 
water zones, that is, when the contrast is small, IP 
chargeability measurements help resolve them. 

16.1.8 Disadvantages 

IP effect being complex in nature, it is difficult to 
interpret the lithology or groundwater condition only 
from IP measurements unambiguously. 

16.1.9 Limitations 

A major limitation in IP is electromagnetic coupling. 
It is an electromagnetic response between grounded 
receiver and transmitter wires used in the survey. The 
effect becomes prominent at large electrode separation. 
Also, the effect is significant for low ground resistivity 
and higher frequencies. It can be mathematically 
removed by observing the phase shift over 3 or more 
frequency decades. A long pulse time combined with 
measurements at long time delay can be used to reduce 
electromagnetic coupling effect. Electromagnetic 
coupling may also be minimized by careful selection 
of electrode array. 

16.2 Mise-s-la-masse 

16.2.1 Purpose 

To map on ground surface the extension of a shallow 
level saturated fracture zone in hard rock which has 
already been encountered in a borehole. 

16.2.2 Principles of Measurement 

In groundwater exploration, when a saturated fracture 
zone in hard rock or a water-filled cavity in limestone 
is indicated by surface geophysical methods, a borehole 
is drilled to tap the fracture zone or cavity. Once the 
borehole has encountered fracture zone or cavity, it 
may be necessary to demarcate its lateral extensions 



for drilling a few more boreholes into the zone of 
interest. To make the operation economic avoiding 
failed boreholes a technique of mapping is used by 
electrically charging the fracture/cavity and tracing the 
developed potential on the ground surface. The 
technique is prevalent in mineral exploration where 
contrast in conductivity of mineralized zone and host 
rock is high. In groundwater exploration, though the 
contrast in conductivity is usually not great, technique 
can be used for lateral demarcation of such fracture 
zones and cavities occurring at shallow level. 

Overall objective is to map developed potential 
distribution on ground surface due to the energization 
of fracture zone/cavity and thus determine its geometry. 
In this case, absolute value of potential is not 
significant. 

16.2.3 Instruments 

Any resistivity meter, preferably with multi-selection 
constant current input can be used to measure the 
potential. Microprocessor based stacking facility for 
rejection of noise is preferred. 

16.2.4 Field Procedures 

In this technique, two types of electrode configurations 
are used. They are the normal (potential) and the 
potential gradient. In normal configuration one of the 
potential electrodes is fixed at a large distance. The 
other potential electrode is moved from one station to 
another along profile lines. In potential gradient 
configuration, the potential difference is measured by 
a potential dipole with a small spacing. Normal 
configuration is preferred as in gradient configuration 
very small potential difference is measured and error 
is added for both the electrodes. Also, speed of 
coverage is better in normal configuration. 

Proper location of infinite or distant current and 
potential electrodes with local geological 
considerations is important. The active current electrode 
is placed against the fracture zone within the borehole. 
Current input should be sufficient for recording a 
considerable potential on the surface. Profile lines can 
be parallel with a spacing of 10 m to 50 m with a station 
interval of 5 m to 10 m surrounding the well. The 
profiles could be radial also, diverging out from the 
well in 8 directions to 12 directions. The pits or holes 
for planting the potential electrodes is made in advance 
and watered to avoid spurious potentials. 

16.2.5 Processing of Data 

If constant current source is not used, potential data 
are normalized to a fixed current input say 1A. The 
potential data can also be converted to apparent 
resistivity data. Data can be taken as profile or contour 
maps of iso -potential or /so-resistivity. 
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16.2.6 Interpretation 

The potential developed depends on thickness of the 
fracture zone, its depth of burial, extent, dip, resistivity 
contrast with the surrounding and presence of 
conductive overburden. In normal configuration 
potential profile broadens with increased depth of 
burial of fracture zone. That is, a deeper fracture zone 
may not show-up on potential distribution map. 
Asymmetry in potential distribution is introduced if 
fracture zone is inclined. Dip of fracture zone is better 
indicated in gradient measurement. Negative peak 
amplitude on up dip side is greater than positive peak 
amplitude on down dip side. Whereas, in normal 
measurements updip direction is indicated by steeper 
slope while downdip side is indicated by gentler slope. 
Presence of conductive overburden reduces the 
magnitude of anomaly while a variation in its 
thickness would complicate the anomaly. Near surface 
inhomogeneities would also affect the anomaly. 
Fracture zone is better defined if it has higher 
conductivity contrast with the surrounding. In contour 
maps elongated contours indicate trend of fracture 
zone. Axis of elongated contour with maximum value 
should be above the fracture zone. 

Location of active current electrode in borehole 
(though it can not be controlled) affects anomaly 
pattern. For a dipping fracture zone if current electrode 
is located in borehole towards its upper part, anomaly 
will be relatively sharp in comparison to that if located 
towards its bottom part. 

In this technique resolution between two fracture 
zones located nearby is poor, that is tracing 
unconnected fracture zone located in proximity of one 
being electrically charged in the borehole is difficult, 
however, better reflected in gradient measurement. 
While interpreting, both resistivity and potential data 
should be studied jointly. 

16.2.7 Advantages 

It can be fruitfully used to ascertain lateral continuity 
of a saturated fracture zone or water filled cavity in 
limestone. The field operation is easy and inexpensive. 

16.2.8 Disadvantages 

A borehole which has encountered fracture zone is 
required to conduct the survey. It cannot be conducted 
in boreholes with pumps running. 



16.2.9 Limitations 

Deeper zones or zones with little contrast in resistivity 
with the surrounding will not be picked-up. Location of 
reference current electrode affects the potential profile 
seriously. Any other fracture zone in vicinity of the 
charged fracture zone would affect potential values. Also, 
near surface inhomogeneities affect the potential values. 

17 REPORTING AND PRESENTATION OF 
RESULTS 

Presentation of geophysical survey results in the form 
of report is very important, as it helps the planners to 
carry out follow-up actions effectively and 
unambiguously. The report should contain summary, 
introduction, objectives, hydrogeological information, 
geophysical approach, quantum of surveys carried out, 
results, discussions, conclusions and recommendations, 
and references. The language should be clear and 
simple without any vague inference. 

The report should invariably contain the 
hydrogeological inferences and the recommendation 
on the objectives or purpose of the survey. It may 
contain technical and physical limitations and 
uncertainties in parameter estimation. Follow-up 
studies or surveys required should be mentioned. 

There has to be a location map at appropriate scale 
showing profile lines, geology and lineaments etc. 
There should be figures showing representative data 
with interpretations. Also, figures to be given as per 
the objective, for example, cross-sections, isopach 
maps, bedrock depth contours, fence diagram, depth 
slicings, correlation with available borehole data etc. 
Maps on aquifer thickness variation, depth to aquifer, 
thickness of protective clay, depth of saline/fresh 
groundwater interface, etc should also be given. The 
interpreted results are generally presented in a tabular 
form. It will not only include the geophysical 
parameters but also the hydrogeological inferences 
in brief and order of priority of the sites to be taken 
up for drilling, with reasons. If the geophysical 
parameters are correlated with hydrogeological 
parameters the graphs as well as map for inferred 
hydrogeological parameters should be prepared. In 
this case regression analysis and confidence limits 
may be shown. Overall, the aim should be to give 
maximum number of utility maps for the objective, 
in a clear and concise manner. 
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